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SYNOPSIS 


The phenomenon of thermophoresis plays a vital role 
in the mass transfer mechanism of several devices involving 
small micron sized particles and large temperature gra- 
dients in the flow field. This dissertation analyses the 
two-dimensional, laminar, thermophoretic flow over an 
inclined plate and over a circular cylinder, and presents 
results for the hydrodynamic, thermal and particle concen- 
tration boundary layers over a wide range of parameters 
for cold, heated and adiabatic wall conditions. Special 
emphasis is placed on the external aerosol deposition 
process . 

The governing boundary layer equations are solved 
using an implicit finite-difference scheme. A marching 
solution procedure is employed whereby the velocity, tem- 
perature and particle concentration distributions at any 
downstream location are determined knowing only the flow 
field upstream of that location. Solution is thus obtained 
in a series of steps covering the complete flow field 


XV 


starting from the stagnation point* While a non-uniform 
grid is used in the marching direction, the grid in the 
transverse direction is self-adapted to the gradient of 
a selected function (velocity or temperature or particle 
concentration) in the flow field. A variable mesh tech- 
nique is employed to preserve the second order accuracy 
of the velocity and temperature profiles in the transverse 
direction . 

For flow resulting from tne impingement of a slot jet 
over an inclined plate, the unsymme trical pressure distri- 
bution on either side of the stagnation point is computed 
by considering the direct impingement of two unequal invi- 
scid jets. The hydrodynamic and thermal boundary layer 
results obtained using the present procedure at zero inci- 
dence are in perfect agreement with the standard similarity 
solution (Blasius flow). Results for local skin friction 
coefficient and Nusselt number in the case of 90° angle of 
inclination of the plate are also in complete agreement 
with earlier results. For angles other than 0 and 90°, 
this is the first study. 

The behaviour of the particle concentration boundary 
layer for the cold plate condition is found to be similar 
to that of the velocity and thermal boundary layers for 
different angles of inclina'tion , of the plate. Near the 
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stagnation point the hydrodynamic, thermal and concentra- 
tion boundary layer thicknesses decrease as p is reduced. 

Away from the stagnation point the situation is reversed 
except for p = 90°. Also, for p = 90°, all boundary layer 
thicknesses are almost constant over half the jet thickness 
on either side of the stagnation point. The presence of 
the stagnation point (at non-zero angles of inclination) 
changes the transverse velocity distribution considerably. 
While it is always positive for p =0, it is negative over 
most of the plate length when p ^ 0. The local skin fri- 
ction coefficient displays a large value close to the 
stagnation point for small values of angle p. For p = 15° 
the local Nusselt number near the stagnation point is 
about ten times that far away from the stagnation point. 
However, this ratio is only about two for P = 90°. 

For cold plate conditions the particle concentration 

at the Wall, C , is found to depend mainly on the wall tempe- 
w 

rature, T_, and the thermophoretic coefficient, K. It remains 
w 

almost independent of the angle of inclination of the plate, 
and the location along the plate as long as K and 
T, , are constant. The wall concentration increases 
as the wall temperature increases or as the thermophoretic 
coefficient decreases; the increase being marginal at 
higher K values. However, with zero angle of inclination 


xvii 


of the plate, the wall concentration drops towards zero 
for thermophoretic coefficients very close to zero. For 
any other angle of inclination the dimensionless wall 
concentration reaches unity at K = 0. For heated plate 
conditions the presence of a critical layer adjacent to 
the wall is demonstrated. It is established that irres- 
pective of the angle of inclination the particle concentra- 
tion in the critical layer near the heated plate approaches 
+ CO when Pr.K ^ 1. The critical layer thickness increases 
as 

(i) the thermophoretic coefficient increases; the 
increase being linear at larger K values 

(ii) the hot Wall temperature increases 

(iii) the distance from the stagnation point increases; 
the increase being appreciable for smaller angles 
of inclination, and 

(iv) the angle of inclination increases. 

The cold wall temperature profiles exhibit a peak within 
the thermal boundary layer at higher Eckert number, Ec , 
values. Accordingly, the particle concentration at the 
inclined plate surface decreases as Ec increases. More- 
over, for smaller inclination angles the dimensionless 
concentration profiles also exhibit a peak value (> l) for 
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Ec > 3. For adiabatic plate conditions the concentra- 

ro 

tion profiles are similar to those for the hot plate 

except that C =0 and the particle concentration gradient 
w 

at the Dlate is zero. 

For flow past a circular cylinder both theoretical 
and experimental pressure distributions are considered 
for the boundary layer analysis. The hydrodynamic points 
of separation obtained in both cases, (105.5° and 81° res- 
pectively), are in perfect agreement with the values 
reported in literature. A comparison of the boundary 
layer results for theoretical and experimental pressure 
distributions at a location $ < 60° from the stagnation 

point shows slightly larger (velocity, thermal and parti- 
cle concentration) boundary layer thicknesses in the latter 
case. However, near $ = 81°, this difference is consi- 
derable. For flow past a cold cylinder also, the wall 
concentration is constant with respect to the nature of 
the pressure distribution over the cylinder (experimental 
or theoretical), and location along the cylinder surface, 
as long as the wall temperature and K are constant. A 
similar observation for flow over an inclined plate 
implies an almost uniform particle 
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concentration 

at any 

two-dimensional 

surface 

in 

laminar flow 

as long 

as T (< 1) and 
w 

K are 

cons tant 

At large Ec 

the peak 

in temperature 

prof ile 

is 


quite appreciable and occurs very close to the cylinder 

surface. Accordingly, the particle concentration within 

the boundary layer increases as high as 20 ^ above that 

of the free stream value (at Ec = 10). With adiabatic 

wall condition (C = 0) this increase is still higher 

w 

(> 40 -/ at Ec = 10). The cold cylinder wall concentration 
reduces as Ec >0 and reaches zero at Ec approaching a 
value of 10. The critical concentration layer thickness 
for a heated cylinder increases as we move away from the 
stagnation point. Eventhough the overall effect of the 
pressure gradient, dP/dX, on the cold wall particle con- 
centration is found to be negligible for Ec = 0, it is 
significant over either a hot wall or a cold wall with 
viscous dissipation. Hence in practical situations 
(where generally Ec 0) it is essential to consider the 
pressure gradient effect. 

From a comparison among the different cases 
studied, smaller angles of inclination when the 
wall is long, and larger angles of inclination 
when it is short are preferred to obtain best results 
on the hot side of heat exchanger equipment where the 


XX 


aim is to minimise particle deposition on the walls. The 
reverse is appropriate for the cold side. Also one should 
be careful about the parameter Pr.K for hot wall applica- 
tions as theoretically C,„ = ■+-« when Pr.K > 1. These con- 
siderations are critical for high temperature gas turbine 
applications as the thermophore tic phenomena could off-set 
remarkably the advantages of transpiration cooling of the 
blade surface. The external small particle deposition 
process provides one of the main mechanisms by which the 
optical fiber preforms are constructed using the Outside 
Vapour Deposition (OVD) process. 
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CHAPTER 1 


INTRODUCTION 


Thermophoresis is the phenomenon by which submicron 
sized particles suspended in a gas acquire a velocity in 
the direction of decreasing temperature. This phenomenon 
causes the aerosol particles to be driven away from a hot 
surface and towards a cold one. This is similar to the 
phenomenon of photophoresis caused by an intense beam of 
light making the particles to move either in the direction 
of radiation flux or against it. The thermophoretic force 
experienced by the small particles is caused by differen- 
tial molecular bombardment giving rise to a radiometric 
force (Green and Lane, 1964), The latter, like other 
interactions between particles and a gas, depends essen- 
tially upon the ratio of the particle radius to the mean 
free path of the gas molecules. The resulting thermopho- 
retic force is proportional to the temperature gradient 
and depends upon the thermal conductivity of aerosol 
particles and carrier gas. Thermophoresis is of consi- 
derable importance for particles as large as 10 pm in 


radius and temperature gradients of the order of 
5 K/mm. 


2 


The phenomenon of therraophoresis is of considerable 
practical importance as detailed by Green and Lane (1964) 
and by Fuchs (1966). Vapour and/or particle mass trans- 
port and/or deposition is central to a variety of engineer- 
ing fields, including materials technology (electronic 
device fabrication, chemical coating of metals, etc.), gas 
clean-up (e.g,, f ilteration) , and the corrosion and/or 
fouling of heat exchangers and gas turbine equipment. In 
the latter case the sharply increasing cost of high-grade 
fuels has stimulated a renewed interest in the combustion 
of coal derived or lower grade residual petroleum derived 
fuels, each associated with inorganic impurities which form 
ash or salt deposits. At the same time, the need for 
higher turbine inlet temperature, required for lower fuel 
consumption^ demands novel, more efficient turbine blade 
cooling schemes, such as film cooling or full coverage 
transpiration cooling. Under these combined circumstances 
thermal diffusion (i.e., thermophoresis for small parti- 
cles) plays an important role in determining mass deposi- 
tion rates, as demonstrated (principally in the absence 
of viscous dissipation and/or transpiration or film cooling) 
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by Goren (1977) and in the more general case by Gokoglu 
and Rosner (1986) o It has been demonstrated that thermo- 
phoresis is the dominant mass transfer mechanism in the 
Modified Chemical Vapour Deposition (MCVD) and Outside 
Vapour Deposition (OVD) processes (used in the fabrication 
of optical fiber preforms) by comparing experimental mea- 
surements and quantitative theoretical predictions 
(Simpkins et al., 1979; Walker et al., 1980; Alam and 
Mehrotra, 1987). Moreover this subject is currently 
of importance in view of its relevance to postulated 
accidents by radioactive particle deposition in nuclear 
reactors. The fact that scrubbing is more effective when 
dusty air is preheated has been ascribed to thermophoresis. 
Another effect of thermophoresis, namely, the repulsion of 
particles from a hot wall and the formation of a dust-free 
layer around hot objects, has also been observed. 

1.1 Literature Purvey 

Movement of particles in a thermal gradient has been 
studied by a number of workers specially in connection with 
the dust-free space surrounding hot bodies. The fundamental 
physical process responsible for the phenomenon of thermo- 
phoresis was first investigated by Maxwell (cf. Kennard, 
1938) , in an attempt to explain the radiometric effect. 
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The extensive literature on thermophoresis contains a 
number of conflicting results, both theoretical and 
experimental. The following gives a brief review of the 
existing important experimental and theoretical studies 
on this problem. 

1.1.1 Thermophoretic Coefficient 

The phenomenon of thermophoresis is caused by the 
so-called radiometric forces, which because of the gaseous 
medium act upon non-unif ormly heated particles present in 
the medium. The radiometric force is developed due to 
the fact that gaseous molecules are repulsed from the 
hotter side of a particle with a greater velocity than 
from the cooler side. That is why particles transmit im- 
pulses directed toward a lower temperature. The velocity 
acquired by the small particles relative to the gas velo- 
city is known as the thermophoretic velocity. The initial 
studies of thermophoretic transport involved the simple 
one-dimensional flovjs for the measurement of this velocity. 
The thermophoretic velocity, w^, is found to be proportional 
to the product of the temperature gradient in the flow 
field and v/t, where v is the kinematic viscosity of the 
gas and t is its absolute temperature. The proportionality 
constant in the expression for is termed as the 



5 


thermophoretic coefficient, K. The value of this 
coefficient depends on the molecular flow regime. A 
measure of this regime is provided by the Knudsen number 
which is the ratio of the molecular mean free path, X , 
to the diameter of the particle, 2Rp. For Rp << X (free 
molecular flow), Waldmann (1959) observed that the thermo- 
phoretic velocity is independent of the particle radius 
and the thermophoretic coefficient has a constant value 
close to 0.75 for elastic collisions of molecules with 
the aerosol particles and 0.54 for collisions resulting in 
condensation and reevaporation. Epstein was the first to 
calculate the thermophoretic velocity in the limit of 
large size particles (R >> X ; slip flow) by solving the 
equations of viscous hydrodynamics (cf. Derjaguin et al., 
1976), With an arbitrary value of thermal slip coeffi- 
cient, = 0.75, the Epstein formula for thermophoretic 
coefficient is 

k 

K = k. 3 (1.1) 

k + k /2 
g P 

where k^ and kp are, respectively, the thermal conductivity 
of the gas and of the particles. Brock (1962) made the 
formula more precise by using an approximate solution of 
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the Navier-Stokes equations for the velocity field near a 
particle while taking into account the particle radius of 
the order of the mean free path, (i.e., R • His 

H 

formula for K is 


K = k 


1 + 


(1+k /2k; + C. :\/R .k /k )(1+2C . >,/R ) 

' P 9 t ' p p m P 


( 1 . 2 ) 


where is the temperature jump coefficient and is the 
coefficient of the gas kinetic isothermal slip. Here again 
k^ is assumed to be 0.75c. The values of range from 
1.875 to 2.48 and values of range from 1.00 to 1.27. A 

review of these values and deductions may be seen in the 
paper by Street (1960). Brock's result, with k^ = 0.75, 
is found not to be in good agreement with the experimental 
data for particles of high thermal conductivity, although 
the discrepancy is much less than that obtained using the 
Epstein result. 


Following a different approach Derjaguin et al. 
(1965 and 1966) reported the following relation for K : 


K = 


IfCt A/Rp.kp/k^ 
l+kp/?kg+C^ k/Rp.kp/kg 


(1.3) 


In the first approximation k^ was assumed to be equal to 
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1.5. Careful experimental measurements with oil mist and 
NaCl aerosol by Derjaguin et al. (1976) indicate that K 
has a maximum value of 0.93 at ~ 1*8* Equation (1.3) 

for K is found to hold good over a large range of ^/Rp 
values and also for a variety of aerosols. Using the experi- 
mental value of thermophoretic velocity as the basis, follow- 
ing values of were derived by Derjaguin et al. (1976) ; 

= 1.13 from Eqn. (1.1), 

= 1.09 from Eqn. (1.2), 
and k^ = 1.07 from Eqn. (1.3). 

Springer (1970) made an attempt to construct an inter- 
polation formula which matches the theory in the entire 
Knudsen number range from near-continuum conditions to free 
molecular (collisionless) limit. However, his formula was 
rather complicated eventhough it agreed reasonably well 
with experimental results. Moreover for the slip flow 
limit, the empirical equation used by Springer employed the 
Maxwell value for k^, which requires for the best fit an 
unrealistic value for C^. A critical review of the data 
and the theoretical expressions for thermophoretic coeffi- 
cients was given by Talbot et al. (1980). An empirical 
expression for K approximately valid for all Knudsen 



8 


numbers has also been given by Talbot et al.^ 

fi j 

K = 2k VigCAl^Bj^ e P 

(1+3C^ x/Rp)(l+2kg/kp+2C^ X/Rp) 


(1.4) 


in which the constants have the values = lo20, = 0*41, 

C, =0.88, k. = 1.17, C = 1.14 and C. =2.18, After a 
series of comparisons with the available data on thermo- 
phoresis , Talbot et al. reported that the empirical rela- 
tion (1.4) aopears reasonably satisfactory for the entire 
range, 0 £ i. Nevertheless, in all the theoretical 

thermophore tic flow analysis reviewed in Section 1.1.2, 
(except that of Gokoglu and Rosner) , it was justified that 
K can be treated as a specified constant of order 1 and 
independent of the temperature in the flow fiela. 

Thermophoresis in laminar flows has been studied 
experimentally by Goldsmith and May (1966), Waldmann and 
Schmitt (1966), Fulford et al. (1971), Derjaguin et al. 
(1976), Talbot et al. (1980), and others. Most of the 
recent experimental determinations of the thermophore tic 
force or velocity have been carried out by one of the two 
methods. One method employs a modified Millikan cell 
wherein a charged particle is held at rest or its velocity 
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measured under the combined action of gravity, thermal 
and electrostatic forces (Waldman and Schmitt, 1966). 

The other method involves the observation of particles 
moving along a narrow channel, either horizontal or verti- 
cal, under the action of a parallel or transverse tempera- 
ture gradient (Derjaguin et alo 1976). A Laser-Doppler 
Velocimeter (LDV) study (by Talbot et al. 1980) of the 
velocity profiles in the laminar boundary layer adjacent 
to a heated flat plate at zero incidence revealed that the 
seed Particles used for LDV measurements were driven away 
from the plate surface by thermophore tic forces, causing a 
particle-free region within the boundary layer of approxi- 
mately one half the boundary layer thickness. For a review 
of the experimental data on therraophoresis Fuchs (1966) 
should be consulted. 

1.1.2 Thermophore tic Flow Investigations 

Thermophoretic deposition of small particles in 
laminar tube flow and in natural convection flow over a 
vertical plate have been analysed by Walker et al. (1979) 
and Epstein et al. (1985), respectively. These two studies 
are similar to that of Goren (1977) which will be elabora- 
ted later in this section. However, the scope of both 
these studies is limited to cold surface conditions only. 
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Assuming a fully developed laminar incompressible flow in 
the circular tube Walker et al. (1979) found that a fra- 
ction of the particles present initially in the stream 
get deposited on the cold tube wall. The effects of 
Brownian diffusion were rigorously included by the use 
of matched asymptotic expansions. However, their similarity 
solution is valid only within a short axial distance of 
the start of constant wall temperature, t^. Using numeri- 
cal methods, the solution was extended downstream and the 
cumulative deposition efficiency at any wall temperature 
and Pr.K value was found to rabidly approach an asymptotic 
limit. Here, Pr is the Prandtl number. The value of 
this limit, E^, is related to Pr.K, the normalised wall 
temperature, [ = t^(t^-t^)], and the wall concentra- 

tion, C , as 
w 

= (Pr.K.C^)/ 0^ . (1.5) 

Both laminar and turbulent natural convection boundary 
layer flow on a cold vertical plate have been considered 
for thermophore tic analysis by Epstein et al. (1985). The 
governing equations reduce to second order ordinary diffe- 
rential equations which are solved numerically for speci- 
fied Pr, K and wall temperature. Results similar to 
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Goren (1977) were obtained for wall concentration in both 
laminar and turbulent natural convection flows. Also it 
was noted that while the particle concentration distribu- 
tion for laminar flow differs notably in shape from that 
for the turbulent flow both profiles converge to practically 
the same oarticle concentration at the wall. Also the 
change in Grashof number, » has no significant effect 

on the wall concentration, C . Epstein et al. derived 

w 

a relation between particle concentration, C and its gra- 

w 

dient at the surface, C'(0), as 


C 


w 




( 1 . 6 ) 


where ©' (o) is the derivative of the normalised tempera- 
ture, 9[=(t-t )/(t -t )], at the wall. Accordingly, it 
was argued that the change in concentration gradient on 
the surface accompanied by a corresponding change in the 
wall temperature gradient keeps the particle concentration 
at the wall almost constant. The simple function 


w 


, _ sloGbt /t 

l-(K.Pr) w «= 

l-(K.Pr)^*^^ 


(1.7) 


was found to relate the computed values of within a 
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relative error of < 3’/ for the values of Pr, Gr , K and 

A 

^v/^oo ir^ves tigated by Epstein et al. Moreover, by re- 
placing the numerical exponent 1.25 by 1,1 in the numera- 
tor and denominator of equation (1.7), the calculated 
results of Goren (1977) for deposition in forced incompre- 
ssible laminar flow over a cold flat plate can be accura- 
tely reproduced. 

Including a volumetric heating term in the energy 
equation the thermoohore tic flow through a circular tube 
has been studied by Morse and co-workers. These studies 
have direct application in the laser modification of 
thermophore tic deposition in MCVD orocess used in the manu- 
facture of optica] fiber preforms. The first work (Morse 
>jnd ClpolLa Jr., lOdJ ) duali< with a tubular 1 low of pai'tl — 
cle laden gas axially heated by the absorption of laser 
energy, neglecting emission and scattering. By a series of 
simplifications they reduced the governing equations for 
temperature and particle concentration into those for a 
one-dimensional problem and solved them for two limiting 
forms of distribution of laser energy throughout the tube 
cross section. Accordingly, it was found that with laser 
radiation of modest power (10 W) , almost all the aerosol is 
evacuated from the central region of the tube causing quite 
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a large aerosol concentration in the vicinity of the tube 
walls A consequence is that in sufficiently long tubes 
the total deposition efficiency can be made to approach 
unity. However, since highly focussed radiation interacts 
with very few particles (causing local heating) which 
through thermophoresis are then repelled from the heating 
zone, the consequent increase in concentration elsewhere is 
negligibleo 

Worse et al. (1985) extended their earlier work with 
the addition of another heat source in the wall with an 
arbitrary temperature distribution apart from the axial 
laser radiation in the circular tube. Solution in this 
case is carried out for a set of property values by a 
finite-difference marching technique. It is found that the 
particle deposition rate increases with the laser intensity. 
The presence of wall heat source further increases the 
deposition efficiency as the heated wall forces the parti- 
cles towards the centre of the tube (region of maximum 
laser intensity) . The buoyancy effect of axial laser 
heating (constant intensity) on the thermophore tic flow 
in a vertical tube is the concern of Wang et al. (1985). 

The particle trajectories are calculated from a one- 
dimensional formulation of the combined natural convection 
problem with constant wall temperature boundary condition. 
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A universal velocity profile has been found that minimizes 
the length over which 99yi of the particles are deposited 
on the wall. However, these results are only qualitative 
since the assumption of no variation of velocity and tem- 
perature in the axial direction is not true. 

A detailed theoretical analysis of laminar thermopho- 
retic flow over a flat plate at zero incidence was pro- 
vided by Goren (1977). A similarity formulation of the 
boundary layer equations was carried out by him in the 
absence of Brownian diffusion whose effect is only to form 
a very thin sublayer even on the boundary layer scale. 

For a cold wall the particle concentration decreases monoto- 
nically from the free stream value as the wall is approached 
and a finite concentration at the wall is attained depending 
upon the wall temperature. With the increase of thermopho- 
retic coefficient, K, the wall concentration is found to 
decrease for the flat plate at zero incidence. For a heated 
wall the presence of a critical layer within the concentra- 
tion boundary layer over a flat plate is proved by Goren. 

He also elaborated the importance of the factor Pr.K in 
determining the concentration profile over a hot wall. For 
zero free stream Mach number (no viscous dissipation), 



15 


the concentration profile becomes identical to the tempera- 
ture profile when Pr.K = 1. Moreover, the particle concen- 
tration at the critical layer near the heated flat plate 
suddenly approaches + «> when Pr.K ^ 1. The thickness of 
the critical concentration layer is found to be a function 
of the wall temperature and it increases parabolically 
with K. The effect of greater than zero (non-zero 
viscous dissipation) has also been studied by Goren for 
the flow past a flat plate at zero incidence. 

The first thermophore tic analysis for an external flow 
in the presence of a non-zero pressure gradient has been 
provided by Homsy et al. (1981). They apply the Blasius 
series solution to the transverse flow past a circular 
cylinder. Results are presented for the particle flux 
coefficients as functions of K and wall temperature for 
Pr =» 0.71. Knowing the characterizing constants of any 
outer potential flow, one can calculate the corresponding 
deposition profile from these particle flux coefficients. 
The local flux profile shown for flow past a circular 
cylinder is parabolic with a maximum value n 0.24 at the 
stagnation point and zero near the cylinder angle, $ = 90^. 
For still larger $ physically impossible negative values 
result. The total deposition rate was found to decrease 
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monotonically with increasing wall temperature. However, 
it should be noted that, in the light of the arguments of 
Epstein et al. (1985) a change in the surface concentration 
flux does not essentially imply a change in the absolute 
wall concentration value as temperature gradient at the 
wall also will be different accordingly (cf. Eqn. (1.6)). 
Alam and Mehrotra (1987) studied the problem of thermo- 
phoretic flow past a circular cylinder specially applied 
to OVD process by using a finite-difference technique. The 
results of this numerical approach compared well with the 
analytical solution of Homsy et al. (1981). 

The common assumptions (to derive the particle con- 
centration equation) in all the foregoing thermophore tic 
flow investigations are : 

(i) the physical properties of the particle-gas mixture 
are those of the gas phase, 

(ii) gas physical property variations are negligible, 

(iii) particles of only a single size are present in the 

gas and the particle concentration is sufficiently 
dilute, 

(iv) the particles are small enough so that in the absence 
of thermophoresis the particles move with the local 
gas velocity. 
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( v) thermophore tic velocity component normal to the plate 
is the only one to be considered, and 

(vi) the effect of Brownian diffusion is restricted only 
to an extremely thin layer next to the surface, and 
can therefore be neglected. 

Quite recently, Gokoglu and Rosner (1984a, b; 1985a, b; 
1986) have brought out a series of theoretical studies 
regarding thermophore tically augmented particle transport 
in free and/or forced convection with variable properties, 
transpiration cooling and/or viscous dissipation. However, 
all their studies are in connection with cold wall boundary 
conditions and without any radiation effects. Forced con- 
vection laminar boundary layers with variable fluid thermo- 
physical properties and particle transport properties are 
treated by them (1986) and mass transfer rate predictions 
are made based on rigorous numerical calculations for solid 
surfaces. While a decrease in mass transfer rate is pre- 
dicted for large sized particles, the flux value increases 
with wall temperature (see also Walker et al., 1979). 
Gokoglu and Rosner (1986) present results only for a flat 
plate and two-dimensional stagnation point flow for laminar 
boundary layers. The mass transfer ratio associated with 
thermophoresis shows a little increase for two-dimensional 
stagnation point flow, although the increase is appreciable 
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(thousand fold) when compared with non-thermophore tic 
flows (especially, for larger sized particles). In 
another work by the same authors (1984a) a rational and 
general, yet simole, correlation suitable for engineering 
predictions and optimizations is developed. It is argued 
that the effect of therraophoresis is to introduce a suc- 
tion - like convective term due to w-j. and a reaction - like 
sink term with an effective rate constant into the mass 
transfer boundary layer equations. The reported laminar 
and turbulent boundary layer correlation in the presence 
of thermophoresis, viscous dissipation and transpiration 
cooling was found in good agreement (within an average 
error of about ly.) with the numerical calculations. Results 
were also presented by Gokoglu and Rosner (1985a) for the 
effect of viscous dissipation on the thermophore tic aerosol 
particle transport across laminar boundary layers. It was 
inferred that as the ratio of cold wall temperature to 
free stream temperature departs from unity, increasing the 
importance of thermophoresis, the dissipation induced modi- 
fication of temperature profiles diminishes, decreasing 
the relative importance of viscous dissipation. It is 
estimated that viscous dissipation should increase particle 
deposition rates by about 2^-/, at t^t^ = 0.8 but only about 
at t^t^ = 0.6 irrespective of particle diameter in the 
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_2 

range 0.3 x 10 ^im ^ £ 0.15 The reduction in 

particle deposition at increased was attributed to the 
thermopnore tic sink operating outside the Brownian diffu- 
sion sublayer. In the case of transpiration cooled sur- 
faces, while real blowing dramatically reduces the rate 
at which particles can diffuse (by Brownian diffusion) to 
the wall, the thermophore tic suction associated with the 
reduced wall temperature offsets some of the fouling rate 
advantages expected due to blowing alone (Gokoglu and 
Rosner, 1984b). It was observed that for a fixed blowing 
rate and wall temperature ratio, the larger particles are 
easily blown away due to their smaller Brownian diffusion 
coefficients. Moreover, larger particles experience a 
stronger counter current in the thinner mass transfer 
Brownian diffusion sublayer, which is comoletely embe- 
dded in the momentum and energy transfer boundary layers. 
It was noted that the effect of increasing blowing rate 
on shifting the particle mass concentration profile (away 
from the wall and reducing the concentration gradient at 
the wall) was appreciable eventhough the corresponding 
profiles of velocity and temperature are not much 
affected for these blowing rates. Also, for the 
larger particles very small blowing rates are suffi- 
cient to reduce deposition rates considerably. 
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but the effect of thermophoresis in offsetting that redu- 
ction is also correspondingly greater. As the blowing 
rate increased, the deposition rate is seen to decrease 
dramatically, but, in principle, no finite blowing rate 
would reduce the deposition rate to zero. A generalised 
'law-of-the wall' representation for the constant property 
turbulent boundary layer structure v>/as obtained with a 
fully developed Couette-flow type approximation and a 
mixing length law modified by blowing or suction (Rosner 
and Fernandez de la Mora, 1982; Gokoglu and Rosner 1985b). 
Since they use law of the wall profile for turbulent 
boundary layer only energy and particle concentration 
equations need to be solved. Turbulent boundary layer 
calculations show that thermophoresis increases the deposi 
tion rate almost in proportion to the temperature ratio, 
the enhancement being dependent on the particle size (for 
large particles the increase is linear with respect to 
t^/ t^) . It was also shown that transpiration cooling 
especially reduces the deposition rate of larger particles 
however, due to the bhermophoresis the reduction is also 
offset more for larger particles and higher blowing rates 
(cooler surfaces). 
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lo2 Present Investigation 

From the foregoing literature survey it is evident 
that thermophore tic analysis for non-zero pressure gradient 
external flows are very limited. The methodology of Homsy 
et al. (1981) has its own limitations because of the use 
of Blasius series for solution. Moreover, they have given 
importance to the surface concentration flux which (as per 
Eqn. (1.6)) is not so relevant a quantity as the practically 
important value of wall particle concentration. The latter 
is true for the study of Gokoglu and Rosner (1986) as well. 
Moreover, Gokoglu and Rosner neglect the streamwise varia- 
tion of all dependent variables, even in the convective 
terms. This simplifies the analysis but is unrealistic. 

In addition, with the exception of Goren (1977) all the 
studies are confined only to the cold wall condition. Non- 
zero pressure gradient in the boundary layer over a body 
occurs commonly in practical thermophore tic applications 
either due to the shape of the body or due to -its orienta- 
tion with respect to the main flow. 

Accordingly, for the present analysis a laminar thermo- 
phoretic flow is considered when an arbitrary pressure dis- 
tribution is imposed on the boundary layer over a two- 
dimensional body. The governing differential equations are 
solved by a relatively simple finite-difference marching 
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technique. The fluid properties are assumed constant. 
However, the computer code can easily handle variable pro- 
perties. Self -adao tive grid generation in the transverse 
direction and prespecified variable mesh size in the march- 
ing direction are used. 

The particle concentration boundary layer adjacent to 
the solid surface is studied for a wide range of parameters. 
Both cold and heated walls are analysed with different 
ratios of the wall to free stream temperature in each- case. 
Moreover, adiabatic wall is also studied. The critical 
layer in the vicinity of a heated wall that has received 
little attention so far in theoretical thermophore tic flow 
studies is analysed in great detail. Assuming the working 
fluid to be air, the Prandtl number is kept constant at 
0.7 for most of the present study. Viscous dissipation 
in the boundary layer is considered with Eckert number 
variation from 0 to 10. The thermophore tic coefficient K 
is treated as a specified constant of 0(1). 

Two specific examples are considered for the two- 
dimensional flow with non-zero pressure gradient in the 
boundary layer. First the flow due to the impingement of 
a slot jet over an arbitrarily inclined plate is analysed. 
The pressure distribution over the inclined plate is 
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computed by considering the direct impingement of two un- 
equal, inviscid jets. Secondly, the flow past a circular 
cylinder is analysed for both theoretical and experimental 
pressure distributions in the boundary layer over the 
cylinder. 



CHAPTER 2 


FORMULATION OF THE PROBLEM 


2,1 Governing Equations 

Consider a gas containing suspended aerosol particles 
exposed to a solid surface as shown in Fig. 2,1. The 
particle concentration is assumed to foe dilute enough for 
the velocity and temperature profiles to pertain to those 
for a particle-free gas. For particles of unit density, 
(i.e., particle density comparable to that of the gas) and 
1 pm radius, this assumption limits the analysis to aerosol 
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concentrations less than about 10 particles per cubic 
centimeter of the gas at normal temperature and pressure 
(Goren, 1977). This is indeed the case for most thermo- 
phoretic applications. The velocity and temperature dis- 
tributions are therefore governed by the usual boundary 
layer equations. 

For a steady, laminar, two-dimensional, incompressible 
flow with constant properties, these equations are 
(Schlichting, 1979) : 
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, . . , du ^ dv ^ 

continuity 57 + ^ = 0 , 


, , du , du<, 

momentum : ^ 


dp d^u 

dy 


dux2 


energy : pCp(u ^ + v ^ + p (^) , 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where (u, v) are the velocity components in (x,y) dire- 
ctions (see Fig. 2.1), p , p and t are, respectively, the 
density, pressure and temperature of the gas, and p, Cp 

and k are respectively, the dynamic viscosity, specific 

y 

heat at constant pressure, and thermal conductivity of 
the gaSo 

The boundary conditions are : 


u(x, 0 ) = 0 , 
v(x, 0 ) = 0 » 

u(x,oo) = Ug(x) , 

u( 0 ,y) = Ug( 0 ) » 

t(x, 0 ) = t (x) for a specified wall temperature , 
w 

or ^ (x, 0 ) = 0 for an adiabatic wall , 
t(x,<») = t^ , 

t(0,y) = t_^ , (2o4) 

where ^^(x) is the potential flow velocity at the body 
surface, t^ is the free stream temperature, and t^ is 
the surface temperature. 



It should be pointed out that the boundary layer 
equations are really not valid at x = 0 In the close 


vicinity of tho stagnation point, the full Navier-Stokes 
and energy equations including the diffusion terms in the 
x-direction must be used. 

Later in Sec. 4.1.3 we will look at the actual limit- 
ing values of X for which the boundary layer analysis for a 
laminar slot jet is valid. Numerically the effect of sing- 
ularity at the stagnation point is confined to a region very 
close to the stagnation point by taking a very small step 
size initially in the X-direction (Hornbeck, 1973). More 
details are given in Sec. 3.3. 

As discussed in Chapter 1, the particles are driven 
away from a hot surface or deposited on a cold one by 
thermophoresis. The velocity acquired by the small parti- 
cles relative to the gas veloci’ty is known as the thermo— 
phoretic velocity. Theoretically the thermophoretic velo- 
city vector, w^, is related to the temperature gradient in 
the flow field by the equation 

= -K ^ vt , (2.5) 

where v = ^ is the kinematic viscosity, and K is the 
thermophoretic coefficient which depends mainly on the 
Knudsen number. A number of formulations are available in 
the literature for the determination of K depending upon 
the regime of flow past the particle (cf. Chapter 1). The- 
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semi-empirical formula proposed by Talbot et al. (1980) can 
be used for a wide range of Knudsen numbers. The numerical 
value of the coefficient K lies between 0o25 to 1.25 in 
most of the cases. For the present work, however, K is 
treated as a specified constant independent of temperature. 
Due to the boundary layer behaviour, the temperature gra- 
dient in the direction normal to the surface, dt/dy, is 
much larger than dt/dx. Also, the gas velocity component u 
is large compared to v. Since the thermophore tic velocity 
components are given by Eqn. (2.5) it is clear that the 
thermophore tic velocity component normal to the surface, 
v^, is the only one of importance. It is given by 

'^t ~ t dy ° (2.6) 

Additional assumptions usually made for thermophore- 
tic analysis (Goren, 1977; Walker et al. , 1979; Talbot et 
al., 1980; Epstein et al., 1985) are : 

(i) Thermal radiation effect is neglected. In fact, in 
the presence of absorbing aerosol, radiation, if 
present, can significantly affect the temperature 
distribution. This is therefore a weak assumption. 

(ii) The particles are assumed to be sufficiently small 

so that the relaxation tiite scale for particle motion 
is several orders of magnitude smaller than that for 
the gas flow. This means that in the absence of 
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thermophores is the particles move with the local 
gas velocity. 

With the foregoing considerations and with the usual 
boundary layer approximations the following conservation 
law is valid for the aerosol particle concentration, c, 
per unit volume of the gas 


ox 


+ 


d[c( v+v^) ] 
__ 


D 


d^c 


(2.7) 


where D is the Brownian diffusion coefficient for the 
particles. This equation is rewritten with the help of 
continuity equation (2.1) as 


u 




^(cv^) 

57 


.2 

D 

dy 


( 2 . 8 ) 


The boundary conditions for this second order diffe- 
rential equation are 


and 


c(0,y) 

= ^'oo» 

c(x,«>) 


c(x,0) 

= 0, 


(2o9) 


if the wall is totally absorbing. Here is the particle 
concentration in the free stream. 
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For micron sized particles in air at normal tempera- 
ture and pressure the Schmidt number (ratio of v to D) is 

5 

of the order of 10 . This means that the Brownian diffu- 
sion sublayer is so thin, even on the boundary layer scale, 
that it does not alter the thermophore tic particle deposi- 
tion rate (Goren, 1977; Walker et al., 1979). Hence the 
only effect of Brownian diffusion is to create an extre- 
mely thin particle concentration sublayer adjacent to the 
wall when the wall is colder than the surrounding gas and 
near a critical layer away from the wall when it is heated. 
Accordingly, the second order diffusion term could be dro- 
pped from Eqn. (2.8), which is rewritten as 


u 


dc 


+ V 


dy 




+ 




dy 


= 0 


( 2 . 10 ) 


It should be noted that by neglecting the Brownian 
diffusion term, Eqn. (2.10) has become first order in y 
and hence only the boundary condition at infinity in (2.9) 
can be satisfied by it. For the cold wall condition, a 
non-zero particle concentration at y = 0 is expected in the 
absence of diffusion. This is determined by solving Eqn, 
(2.10) using only the boundary condition at infinity in 
(2.9). Essentially, this concentration is the one at the 
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outer edge of the very thin Brownian sublayer adjacent to 
the solid surface. 

For a heated wall, the presence of a critical layer 
near the wall in the absence of Bro'jvnian diffusion can 
be ascertained from a similarity analysis (Goren, 1977) » 
While Goren' s analysis is for a flat plate at zero inci- 
dence, we carry it out for thermophoretic flow over a 
wedge . 

The similarity variable, T), for flow past a wedge is 
taken to be 


n 


y 




( 2 « 11 ) 


where m is the wedge constant related the wedge angle Tip, 
by 


~ 2 ^ * ( 2 . 12 ) 

In the neighbourhood of the stagnation point, the poten- 
tial flow velocity Ug(x) over the wedge is related to the 
length coordinate, x, measured from the stagnation point 
by 

U-(x) «: x”^ . 

S ' 


(2.13) 
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For the similarity solution, continuity equation is 
taken care of by defining a stream function 'i['(x,y) as 

f(’l) • (2-14) 

This leads to velocity components u and v as 


u = Ugf (t)) , 


(2„15) 


V = 



vu 
s 

X 


[f + 


m—l 

m+1 




f'3 , 


(2.16) 


and the momentum equation reduces to the famous Falkner- 
Skan equation 

f". ^ t f,. + ^ ^ Q ^ 

with the boundary conditions 

f(0) = f*(0) = 0, f ‘ («) = 1 . (2.18) 

Here primes denote differentiation with respect to ri. 

For the thermal boundary layer over an isothermal 
wall, it is known (Schlichting, 1979) that the similarity 
solution is possible only for flow over a flat plate if 
viscous dissipation is included. However, if viscous 
dissipation is neglected, similarity solution exists for 
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flow over a wedge even when the wall temperature varies 

3 S 


[t (x) - t ] oc X 


n 


(2.19) 


For an isothermal wall, of course, n = 0. In this case 
the energy equation reduces to 


= -Pr f(T)) , 


( 2 , 20 ) 


where t' = dt/dr;. Thus the temperature is a function of 
the similarity variable t). With this, the thermophore tic 
velocity is given by 


^t ^ 


m+1 „ t * 

X ^ t ’ 


( 2 . 21 ) 


from which it follows that 


^^t m+1 ^s t‘^ „ t",, 


( 2 . 22 ) 


The simplification of the particle concentration equation 
(Eqn. ( 20 IO)) using the foregoing relations yields 


(K 1’+ f) c’ - (K — - K.Pr f)] c = 0 , (2.23) 



33 


where c' = dc/dt). This first order differential equation 
in c can be integrated with free stream boundary condition 
for the particle concentration c{^) = c to yield 

CO 


£- 

c 


= exp [• 


/ 

n 


t' / j, t* 

" t — 


+ K Pr f) 


(K + f) 


dri] 


(2.24) 


The equation for particle concentration of Goren (1977) 
without the viscous dissipation term is exactly the same as 
Eqn. (2.24) above. Accordingly, the arguments of Goren 
about the behaviour of the concentration boundary layer 
over a flat plate are applicable to the flow over a wedge 
as well. Thus, a critical layer, exists near a heated 
wedge as well. There is no reason to believe that such 
would not be the case for flow over any heated surface. 

There is a value t), say , at which the denominator 
of thr integrand in (2.24) goes tn voro and at which a 
singularity in the dust concentration profile in the 
absence of Brownian diffusion can be expected. The dust 
concentration near the singularity depends ufjon whether 
Pr.K is less than or greater than unity, for the numerator 
in (2.24) also changes sign at some value of ri. For 
Pr.K < 1 the numerator changes sign at a value of g larger 
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than , and the dust concentration very rapidly approa- 
ches zero at the critical layers For Pr^K > 1 the numera- 
tor changes sign at a value of -g smaller than ti , and the 

c 

dust concentration rapidly becomes infinite. Thus depend- 
ing upon Pr.K. ^ l»c •* + “’3sy*»0 for a hot wall. 

2o2 Strategy for Solution 

Analytical solutions to Eqns (2.1) to (2.3) are rare 
due to the presence of nonlinear terms. As is well known, 
an analytical (similarity) solution for Eqns. (2,1) and 
(2.2) with appropriate boundary conditions in Eqn. (2,4) 
is possible but for a very restricted class of non-zero 
pressure gradient, dp/dx. The familiar examples include the 
stagnation point flow, the converging channel flow, and the 
wedge flows, and are well documented in almost every text 
on Fluid Mechanics. Even these solutions, however, are 
valid only in a narrow region close to the stagnation point 
due to the limitations of the function representing the 
pressure gradient term. To obtain the similarity solution 
for the energy equation (2.3) as well, additional constraints 

are required on the variation of t^. Thus an analytical 
solution for a general dp/dx is not possible. 
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The system of partial differential equations (2.1) 
and (2.2) is of parabolic type which can be solved easily 
by finite-difference techniaues while marching downstream 
in the x-direction. Moreover, the nonlinear inertia terms 
in the momentum equation present no particular problem to 
the finite-difference method. Since the energy equation 
is decoupled from the continuity and momentum equations, 
it can be solved once u and v are known. Finally the 
particle concentration equation (2.10) can be solved. 

2.3 Dimensionless Representation 

Before obtaining a numerical solution the governing 
equations are placed in dimensionless form by using the 
following dimensionless variables : 

X = x/L 

Y = (y/L) Re^/^, 

U = u/u^ , 

V = (v/uj Re^"^^ , 

T = t/t^ , 

P = (P-Pj/( PU^) r 

= v^(lA) Re“^/^ , 

C = c/c 


and 


(2.25) 
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Here L is a characteristic length, and Re is the 
Reynolds number defined as Re = u^ 'L/v , For the incli- 
ned plate considered here, L is the width of the impinging 
jet (Fig. 2,2), but for the circular cylinder, L is the 
radius of the cylinder. In addition the following dimen- 
sionless numbers are introduced to further simplify the 
energy equation : 

Prandtl number : Pr = p C /k 

Eckert number : Ec = u^/C t . (2.26) 

oo p oo 


All these quantities are now substituted in the 
governing equations (2.1) to (2.3), (2.6) and (2.10) 
resulting non-dimensional equations are 


c ontinui ty 


momentum 


energy 


• ^ dY ^ 


. U ^ + V 
* ^ dX ^ 


d^U dP 


^ “ dX » 

tf ^ dT 1 d^T ^ 

= ^ ^ ^ 3y = P? ^ » 


thermophore tic 
velocity 


. V - - 
. Vj - T Fy * 


particle 

concentration 


.p dCC V ) 

* ^ ^ 


. The 


(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 
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The general boundary conditions (2,4) and (2,9) 
rewritten for the dimensionless problem become : 

U(X,0) = 0 

V(X.O) = 0 

U(X,=o) = Ug(X) , 

U(0,Y) = Ug(0) , 

T(X,0) = (isothermal wall), 

or =0 (adiabatic wall) , 

T(X,«) = 1 , 

T(0,Y) = 1 , 

C(X,~) = 1 , 

and C(0,Y) = 1 » (2,32) 

where = u„/u 

s s » 

2,4 Pressure Distribution 

The estimation of the correct pressure gradient in 
the boundary layer is essential before solving the momen~ 
turn equation. The pressure gradient is usually computed 
by considering inviscid flow over the body surface. Two 
examples are considered for analysis here. The first one 
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deals with the oblique impingement of a slot jet on a 
plane wall. The oressure distribution in this case is 
unsymme trical about the stagnation point. The second 
example is the flow/ past a circular cylinder which has a 
symmetrical pressure distribution with respect to the 
stagnation point. 

2o4.1 Flow Over an Inclined Plate 


In order to estimate the pressure gradient in this 
case, a direct impact of two unequal jets is considered 
as shown in Fig. 2o2o It is assum.ed that the flow is 
inviscid and the outgoing jets (B and D) are of equal width 
L. In this general case, the equation relating the com- 
plex variable z with the complex velocity | is (Milne- 
Thomson, 1960) : 


u 


z = 


% 


h, t f 

[— ln(l - 4“) + In (1 - —-) 




( 2 . 33 ) 


where , a^^ = , 

e”^^ , 

b2 = e^^ , 


and 


(2„34) 
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Also the following relations are applicable between the 
widths of the jets. 

h^ = L (1 - cos , (2.35) 

h 2 = L (1 + cos p) . (2.36) 

Substituting Eqns. (2.34) to (2.36) into Eqn. (2.33) and 
simplifying results in the following relation : 

TtZ + ln( 1 + £) -In (1- c) - cos p In (1- G^) 

+ ln(l - £ e^^) + In ( 1 - G e""^^) = 0, (2.37) 

where Z = z/L and 0 = |/u^ . 

Without loss of generality, the configuration shovm in 
Fig. 2.2 can be regarded as the jet from D impinging on the 
plane wall along the line AC, which is at an angle p to the 
oncoming slot jet (Fig. 2.3). Hence the dimensionless com- 
plex velocity, £ , obtained by solving Eqno-(2.37) can 

be used directly in order to compute the pressure distri- 
bution over the plane wall AC. In fact, since we are 
interested only in the potential flow velocity Ug(X) at the 
plate, we set Z = -X and £ = U in Eqn. (2.37) to obtain: 

o 
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1-U^ 1+U^-2U cos P 

TiX = -In (T-rr) + cos p In ( § §-,5 ) 

l+Us 1 - 

s 

U sin p 

+ 2 sin p tan (-^ ) . (2.38) 

1-U, cos p 
s 

Thus for a given value of angle p, U (X) is a solution 
of Eqn. (2.38). This equation is identical to Eqn. (3) 
of Miyazaki and Silberman (1972') in the limit H 
The pressure distribution and the oressure gradient in 
the boundary layer over the plane wall are given by 


2P 



(2.39) 


and 



dU 
s 

dX 


(2.40) 


This pressure gradient is used in the momentum equation 
(2.28). 

The pressure distribution in the boundary layer at 
various angles of inclination of the date is given in 
Fig. 2,4. The curves are shown at 15° intervals from 
0 to 90°. As observed experimentally by Garden and 
Akfirat (1966) in Fig. 7 of their paper, the pressure dis- 
tribution on the plate for normal impingement (p = 90°) 
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is symmetrical about the stagnation point (X=0) . It is 
clear from Fig. 2.4 that except for the 90° angle the 
distribution is unsymmetrical about the stagnation point. 
Also the symmetrical distribution for p = 90° cuts 
across the unsymmetrical distribution for smaller angles 
at some positive values of X. The presence of the stag- 
nation point for angles other than zero considerably 
modifies the distribution especially near the stagnation 
point. In the positive X-direction, pressure reaches a 
nearly zero value at X ^^ 3.0 for all angles, but in the 
negative X-direction, zero pressure is obtained at 
different X-locations depending upon the value of p . The 
pressure profiles in the negative X-direction come very 
close to the P-axis for smaller angles (P < 30°) as seen 
from Fig. 2.4. Thus for 6 = 15°, P « 0 for X < -0.3. 
Moreover, for all values of p, P 2:0 forX^ 3.0. The 
reason for this lies in the smaller thickness (h^^) of 
the outgoing stream on the (X < 0) side of the plate at 
smaller values of p (see Eqn. (2.35) and Fig. 2.3). In 
the limiting case of p -* 0, ■* 0 and we get the flow 

over a flat plate at zero incidence. The pressure gra- 
dient with respect to X is shown in Fig. 2, 5. Note the 
sharply changing values of dP/dX near the stagnation point 
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for small angles of impingement of the slot jet. 

Figure 2.6 illustrates the potential flov^ velocity at 
the plate as a function of distance along the plate with 
P as a parameter. The effect of pressure distribution 
on the plate is very clearly reflected in these curves. 

2.4.2 Flow Past a Circular Cylinder 

This configuration is shown in Fig. 2,7. The invis- 
cid flow analysis (superimposing a uniform stream with a 
doublet) yields the following expression for non-dimen- 
sional velocity over the cylinder 


U =2 sin X . 
s 


(2.41) 


Here X corresponds to the location on the cylinder sur- 
face from the stagnation point ( # = x/L = X) . The 
pressure gradient follows from Eqn. (2.41) as 


dX 


-U 


dU 
s d)r 


» 


or ^ = -4 sin X cos X . (2.42) 

Numerical solution of the boundary layer equations 
(2.27) and (2.28) for incompressible flow past a circular 
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cylinder with the pressure gradient calculated using 
Eqn. (2.42) predicts the separation around $ = 105° 
(White, 1974). However, the broad wake caused by bluff- 
body separation is so different from the potential flow 
that it alters greatly everywhere. The experimental 
results for the flow past a circular cylinder fit the 
polynomial (White, 1974) 

Ug = 1.814 X - 0.271 - 0.0471 X^ , (2.43) 

which is quite different from the potential flow. The 
expression for pressure gradient derived from Eqn. (2,43) 
is 

II = -(1.814 X - 0.271 X^ - 0.0471 X^)* 

(1.814 - 0.813 X^ - 0o2355 x"^) . (2.44) 

The finite-difference solution of Eqns, (2.27) and (2.28) 
by Smith and Clutter (1963) with dP/dX given by Eqn, (2.44) 
places separation at a realistic f = 80°, The present 
analysis is carried out with both theoretical and experi- 
mental pressure distributions (Eqns, (2.42) and (2.44) 
respectively). A comparison of the results obtained is 
given in Chapter 5, 



The theoretical and experimental pressure distribu- 
tions in the boundary layer at different locations on 
the circular cylinder are shown in Fig. 2.8. From this 
figure it is clear that the pressure gradient, dP/dX, 
changes sign at 90° and 70°, respectively, for the theore- 
tical and experimental distributions. Hence a point of 
inflection is expected in the laminar velocity profiles 
at angles larger than theseo 




Fig. 2.1 Typical two-dimensional boundary layer configuration 
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Potential flow velocity at the plate 






CHAPTER 3 


FINITE-DIFFERENCE SOLUTION 

3.1 Discretization of the Governing Equations 

In order to represent the dimensionless differential 
equations (2.27) to (2.31) in finite-difference form, a 
rectangular mesh is superimposed on the flow field. Indi- 
ces (j, k) are used to indicate position in the (X, Y) 
directions respectively. The finite-difference coordina- 
tes are chosen to correspond to the spatial coordinates 
(shown in Figs. 2.3 and 2.7) so that j = 0 corresponds to 
X = 0 and k = 0 corresponds to Y = 0. Accordingly, (0, 0) 
represents the stagnation point both in terms of indices 
and spatial coordinates. A small change in the X-coordinate 
(AX) increases j by 1 and similarly an increment by AY in 
the Y— coordinate increases k by 1. The wall is represented 
by k = 0 and edge of the boundary layer by k = n+.l. 

The discretization of the governing equations is made 
in such a way that the solution can be carried out by a 
marching procedure in the X— direction. The momentum, energy 
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(shown in Figs. 2,3 and 2.7) so that j = 0 corresponds to 
X = 0 and k = 0 corresponds to Y = 0. Accordingly, (0, 0) 
represents the stagnation point both in terms of indices 
and spatial coordinates. A small change in the X-coordinate 
( AX) increases j by 1 and similarly an increment by AY in 
the Y-coordinate increases k by 1. The wall is represented 
by k = 0 and edge of the boundary layer by k = n+1. 

The discretization of the governing equations is made 
in such a way that the solution can be carried out by a 
marching procedure in the X-direction. The momentum, energy 


54 


and concentration equations are solved implicitly at 
every X-location, whereas the continuity equation is 
solved in a stepwise manner. Following Hornbeck (1973) 
forward differencing in the marching X—direction and 
central differencing in the Y-direction are employed in 
Eqns. (2.28) and (2.29), whereas only forward differen- 
cing is employed in Eqns. (2.27) and (2.30). The finite- 
difference form of these equations is written as : 

Continuity Equation : 

tlrktl ~ » k+1 "^j+lrk+l ^i4-l,k _ Q ^ (3.1 

AX AY ' * 


Momentum Equation : 


U. , U. , U. , , - U. , , , 

U l + A sX Ir-k + V 

^j+l,k AX j+l,k 2( AY) 


^1+1, k+1 “ ^^.i+l,k ^ ^1+1, k-1 _ ^ 
( A Y)^ 


(3.2) 


Energy Equation : 


n ’^.i + l»k ~ ^j ,k y .Il+lt.K.7,3-. 

^j + l,k AX j+l,i< 2( AY) 

- !_ ^1+l,k-fl ~ ~^i+l,k ^ '^1+1, k-1 


Pr 


( A Y)' 


U 


- U. 


+ Ec [ " . 1 t 2 


f 


(3.3) 
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Thermophore tic Velocity : 


M = _ L- r ^1+l.lc+l “ '^1+1. k . 

Vj AY T ■' 


(3.4) 


These discretizations are second order accurate in aY 
for U and T. 

While obtaining the solution by the Thomas algorithm 
(Appendix A), the sufficient conditions to be satisfied 
by the tridiagonal coefficient matrix elements are 
(Roache, 1982) : 


( i) the diagonal elements should be positive, 

(ii) the off-diagonal elements should be negative or zero, 

(iii) the magnitude of the diagonal element should be 

greater than or equal to the total magnitude of the 
two off-diagonal elements in each row. 

In order to satisfy these conditions, the particle concen- 
tration equation (2.31) is discretized in four different 
ways depending upon the signs of V and V-p, and the wall 
temperature condition as follows : 


Cold Wall (V < 0) : 


^j+l»k 


S+l,k+l "* S+l,k 




— ilis + V. , . 

AX J+l,k 

Vt _ r 

^i+l.k+1 ^1+1, k+1 j-hl,k i±Us = 0 , 


(3,5) 
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Cold Wall (V > 0) : 



Hot Wall (V < 0) r 



Hot Wall (V > 0) : 



These discretizations are first order accurate in' AX and 
AY for C. It should be mentioned here that V does not 
change sign at all Y-locations at a particular X. Values 
of V close to the plate become positive first. Hence it 
is difficult to strictly satisfy the three conditions 
mentioned above. However, satisfactory results are 
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obtained for all the cases studied here when the sign of 
V at 25 AY below Y = <» is used as the criteria for chang- 
ing the scheme for concentration calculation (Eqns. (3.5) 
to (3.8)). 

The finite-difference forms of the boundary condi- 
tions (2.32) are obvious except for the adiabatic wall 
condition. A difference form for the gradient which is 
consistent in truncation error with the discretised energy 
equation (3.3) is 

aT(x,o 
57 

where T.,, is the unknown surface temperature. 

j+1,0 

3.2 Method of Solution 

Equations (3.1) to (3.3) and (3.5) to (3.8) are re- 
written in forms which are more readily adaptable for 
soluti on . 

Continuity Equation : 

r AY 

''^j+l,k+l " ^j+l,k ^^j,k4-l '^j+ljk+l^ aX 


■^"^1+1,0 ^ ^"^i+l,! 


i+1,2 


(3.9) 


(3ol0) 
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Momentum Equation 


V. 


[- 2^*^) ■ (ay) 2^ “j+l,lc-l + f 


u 




i 1 T T 

( A Y)^ j+1 ,k+l 


j+l J , 

2 

— . 1 IT 

AX + 

( A Y) 

2 J u 


1 I 

dP 

AX 

'J » , 

J »k 

■ dX 


(3.11) 


Energy Equation : 


r _ _ i 1 T . r 1 + 1 „ k. 2 T _ 

2t A Y) ^j+l,k-l L AX + . ,,,2^ 


Pr( A Y) 


2-' j+l,k 


1 — ^ ^j+l>k p 

Pr( A Y)^ j+l>k+l AX j »k 


‘-Hay) 


(3.12) 


Concentration Equation : 


Hot Wall (V > 0) : 

V ^ 

r .l+lxK _ 

*- AY 


T. 


AY J ^ 


ji-l,k~l 


^ AX ^ AY ^ AY J j+l,k 


AY 


AY 


r 

AX ^ j , k 


(3.13) 
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Hot Wall (V < 0) : 


1+1 ; k — 1 _ 

AY 3+l,k-l 


[ 


U 




AX 


V 

1 , K ^ .1+1 T k-| Q 

AY AY -* j+l,k 


V. - , 

, 1+1 » k p 

AY j+l,k+l 


1 'i' 1 # k 


Cold Wall (V > 0) : 


(3.14) 


V 

i+l ♦ k p 
AY j + l,k-l 


+ 


itii 
A 


+ 


V 

llthl 

AY 


T 

1+1 .k -i p 

AY 3+1, k 


1+ltk+l r 
AY 3+l,k+l 


^■3+1 e 

AX 3,k 


(3.15) 


Cold Wall (V < 0) : 


r !i t . X + 

'-AY AX 


V 


T., 


iiiiiS] c 

AY •' 3 + 1, k 


V T IJ 

+ + j+i>k+i i p ^ _ 

'-AY ^ AY J 3 + 1, k+1 AX 


i+1 1 k p 


(3.16) 


As mentioned earlier, the solution for the variables 
U, V, T and C is carried out in a marching procedure start- 
ing from the stagnation point (X=C). The discretized momen- 
tum equation (3.11) written for k = 1 (1) n leads to a set 
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of n nonlinear equations for n unknown values of U at 
the X-location (j+1). This set is reduced to a linear 
set by replacing the coefficients U. , , and V. ^ , in 
(3*2) by their known values at the orevious iteration 
(indicated by the superscript* in £qn. (3.11)). The 
linearised tridiagonal set is solved iteratively using 
the discretized continuity equation (3.10) to update V 
at the location (j+1). 

Once values U and V have been determined at (j+1), 
the discretized energy equation (3.12) written for 
k = 1 ( 1) n leads to a tridiagonal set of n linear equa- 
tions that can be easily solved for n unknown values of T 
at the location (j+1). If the adiabatic wall boundary 
condition is used the additional unknovm wall temperature 
requires an additional equation for solution which is 
obtained from Eqn. (3.9). In this case, the last element 
in the first row of the coefficient matrix requires to be 
eliminated to preserve the tridiagonal character of the 
system of equations. This is done using the Gaussian 
e limination . 

The temperature values computed at the (j+1) location 
are substituted in Eqn. (3.4) to determine the thermophore- 
tic velocity at this location. 
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Solution can then be obtained for the discretized 
concentration equations (3.13) to (3.16) by substituting 
the calculated U, V and Vj values at the (j+l) location. 

For a hot wall, Eqn. (3.13) or (3.14) is selected to repre- 
sent the conservation of the particles in the flow field 
depending, respectively, upon whether V is positive or 
negative. Equation (3.13) written for k = 1 (1) n leads 
to a bidiagonal set of n linear equations whereas Eqn. (3.14) 
leads to a tridiagonal set of the same number of equations. 

The n unknown values of C at the location (j+l) are obtained 
by solving this set of n equations. 

There is a critical layer present near a hot wall. As 
discussed in Sec. 2.1, the particle concentration gradient at 

the edge of this critical layer is very large. Thus, C > + 

> 1 For 

as Y > 0 depending upon the value of PrK being < 

gases, of course, PrK > 1 is of academic interest only since 

Pr 2 ; 0.7 and K < 1. Since the thickness of this critical 

layer is unknown apriori, the large concentration gradient at 
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the edge of this layer can be satisfied by imposing C = + co 
at the wall depending upon the value of PrK. For PrK < 1, the 
physical solution corresponds to C = 0 within the (particle- 
free) layer, and the particle concentration rising sharply at 
the edge of this layer. Unfortunately, an asymptotic solution 
of Eqn. (2.31) lends no help since this equation reduces to 

dc/dy > 0 as Y > «>; a well known fact. At any other 

value of Y, simplification of Eqn. (2.31) is not possible. We 
checked our solution in several ways so as to ensure that 
we had a unique solution. For a flat plate at zero incidence, 
it matched perfectly with Goren’ s results (Goren, 1977). More 
details are available in Sec. 4.3.3. 

For the cold wall condition, the particle concentra- 
tion is finite at the wall as observed earlier (see 
Sec. 2.1). In the case of negative V, Eqn. (3.16) is 
selected as the discretized concentration equation. 

This equation v>/ritten for k = 0 (1) n leads to a bidiago— 
nal set of (n+l) linear equations for the (n+l) unknown 
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values of C at the X~location (j+1). In the case of 
positive V, Eqn. (3.15) is selected to represent the 
particle conservation. However, this equation can not be 
written at k = 0, where again Eqn. (3.16) has to be wri- 
tten and multiplied by a negative sign throughout in 
order to satisfy the conditions given in Sec. 3.1. This 
equation and Eqn. (3.15) written for k = 1 (1) n together 
give a tridiagonal set of n+1 linear equations for the 
n+1 unknown values of C at the (j+1) location. 

Thus the U, V, T, and C values are obtained in 
sequence at every (j+1) location during the marching pro- 
cedure along the X-direction. 

3«3 Computational Details 

The computer programme takes care of the fact that 
the boundary layer thickness has to be increased as we 
move away from the stagnation point. Thus the value of n 
is chosen so as to ensure that there are atleast 3 to 4 
points for k £ n for which U for velocity boundary 

layer. Similar adjustments are done separately for tem- 
perature and concentration boundary layers in which case : 

I 

the comparison is done with T and C respectively at the 
edge of the boundary layer. 
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For computational efficiency, but without sacrificing 
accuracy, the step size AX is changed systematically as 
we march along the X-direction. Initially while starting 
from the stagnation point AX is very small (of the order 
of 10 ) . The AX value is approximately doubled at an 

interval of about 10 marching steps. However, beyond 
jx| = 0.1, marching is done for a large number of steps 
with the same AX (as AX has become sufficiently large). 
Thus in the span of X from 0 (stagnation point) to 3 
(far away from the stagnation point), eleven different 
step sizes (A X) were used for computation. The set of 
AX and the number of steps marched with each AX are 
given in Taoie 3,1, For a snail enough value of the 

A ***5 

initial AX, such as 10 used here, it Is found that 
the sif.gularity at tne stagnation poin L is confined to 
aoout 5 (Ax) on eith„-r side of xhe stagnation point. 

While AX can be assigned any value without difficulty, 
it is easier to keep AY uniform due to central differen- 
cing employed for several terms in the Y-direction. How- 
ever, due to large gradients of velocity, temperature and 
concentration (hot wall) around the critical layer, AY 
must be kept quite small near these regions. Then if AY 
is uniform across the boundary layer thickness, the number 
of simultaneous equations increases excessively. Solution 
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Table 3.1 Marching Pattern along X-direction 


Serial 

No. 

jAXl 

Number of 

steps 

marched 

Final 
jxl value 

1 

0.00001 

20 

0.0002 

2 

0.00002 

10 

0.0004 

3 

0.00005 

12 

0.0010 

4 

0.00010 

10 

0.0020 

5 

0.00020 

10 

0.0040 

6 

0.00050 

12 

0.0100 

7 

0.00100 

10 

0.0200 

8 

0.00200 

10 

0.0400 

9 

0.00500 

12 

0.1000 

10 

0.01000 

90 

1.0000 

11 

0.02000 

100 

3.0000 
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of such large number of equations not only requires exce- 
ssive computer time but also involves large round-off 

error. A practical solution is to use a fine mesh size 
(ay) in regions of large gradients and a relatively 

coarse grid away from it. This requires modification of 
the discretized equations (as explained in Appendix B) 
at locations where the mesh size changes. The sudden 
change in the concentration gradient especially near the 
critical layer (in the case of hot wall) makes the assign- 
ment of different AY in different regions unpredictable. 
This situation is averted by the use of an efficient self- 
adaptive grid scheme along the Y-direction (Section 3o4). 
Also this makes the programme more efficient apart from 
yielding more accurate results at the expense of a margi- 
nal increase in the computational time. 

3.4 Adaptive Grid Generation 

The self-adaptive grid scheme employed is essentially 
the same as that suggested by Nakahashi and Deiwert (1984, 
1986, 1987). Based on the variational principle, a spring 
analogy is used in this method to redistribute the grid 
points in an optimal sense in order to reduce the overall 
solution error. Only one user specified parameter (for 
one-dimensional adaptation) denoting the ratio of maximum 
to minimum grid spacing desired is used to define all the 
important constants, thereby minimizing empiricism and 
making the method self-adaptive. 
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In the one-dimensional adaptation procedure, grid 
points in one direction, say, the Y— direction are re- 
arranged so as to equidis tribute some positive weight 
function w(Y) over the field, i.e., 


'^k+l 

/ w(Y) 


Y. 


dY = constant . 


(3.17) 


In discrete form this equation is represented as 


AYk Wk = constant , (3.18) 

where AYk = Yk+1 “ \ over the interval 

AYk* For the equidistribution scheme with n number of 
grid intervals over a length £ in the physical space, 
grid spacing is given by 

AYk * (3*19) 

As suggested by Nakahashi and Deiwert (I'PS?), the 
expression for weight function, w is chosen to be 

Wk~l+Afk» (3. 20 ) 

where fk - ^ ^ic""^niin^ ^max ^min^ 
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Equation (3.20) has been normalised so that the first 
constant is unity. When the constant A is zero, grid 
points are equally soaced, is any non-negative fun- 

ction (at the grid point k) based on the solution of the 
problem. It is usually taken to be the gradient of the 
flow solution at the grid point k, A and B are positive 
constants, and f . and f„,^ are the minimum and maximum 
values of the function f. 

The constant A controls the ratio of maximum to 
minimum grid spacing desired. This ratio 
Max jAYj^l/Min is equal to the ratio of the weight 

functions, Max | 1 /Min j | , as evident from Eqn. (3.19). 

Therefore, 

AY 

A = - 1 , (3.21) 

^rain 

Thus (A+1) is initially specified as the desirable ratio 
of maximum to minimum grid spacing. The constant B is 
chosen such that the minimum grid spacing MinjAYj^l is 
equal to the specified AY . , and this can be achieved by 

finding solution to the equation 

F(B) = MinjAYj^l - AY^^^ = 0 . (3.22) 
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The Newton-Raphs on method can be used to solve this equa- 
tion iteratively using the following equation 


bN+1 





(3.23) 


where 


= 


F(B^) 

dF(B^) 

3b 


MinlAYj^i - 

^MinlAYj^i 

Tb 


ay . 

min 


(3.24) 


The derivative in the denominator of Eqn. (3.24) can be 
evaluated using Eqns. (3.19) and (3.20) as 


dMin jA Yj^ j 
^B 


n 

-A 2 
i=l 


[AY^(6 


ki 





_B _ 

) f^ In f^] , 


(3.25) 


where is the Kronecker delta, and and corres- 

pond to MinjAYj^j. The iteration for B starts with an 
initial guess (B*^) and the Eqns. (3.19) through (3.25) 
are solved to find a new value of B, namely B^. The 
solution function, the gradient of which is used in 
Eqn. (3.20) as f, is interpolated on the new grid result- 
ing from Eqn, (3.19). Lagrangian interpolation is used 

for this purpose. With the new value of B, the whole 
process is repeated till a reasonable low value of jABj 

in Eqn. (3.23) is achieved. 

A FORTRAN code is written for the self-adaptation of 
grid in the Y-direction following the above steps. This 
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subroutine is called at pre-specif ied X-locations in the 
marching direction so as to readapt the grid points along 
the Y-direction. Any one of the boundary layer solutions 
(U, T or C) can be selected to provide values of f as the 
Y-gradient of U, T or C for adaptation of the grid. How- 
ever, for the hot wall case the particle concentration C, 
which has a rapid variation near the critical layer thick- 
ness, is more appropriate. Section 4.3.3 provides a token 
comparison of concentration profiles with and without grid 
adaptation . 

The minimum value of AY is set at 0.00b and the 
maximum at 0,05. The self-adaptive grid generation tech- 
nique distributes values of AY within this rangeu As many 
as 500 steps in the Y-direction aie required for some val- 
ues of angle p at the farthest downstream location X = 3, 
This discretization is used following considerable experi- 
mentation keeping the error incurred and reasonable computer 
time in view. The test case for error analysis is the 
Blasius solution for the fiat plate at zero incidence. Our 

values of U and V for ^ = 0 are found to match exactly v;ith 

^ —4 

the Blasius solution for all X > 10 , Similar results are 

obtained for 0 for v^hich the exact solution corresponding 
to p = 0 is given by Schlichting (1979). During the test- 
ing phase, values of local skin friction coefficient and 
Nusselt number for normal iopingetment of the slot jet on the 

plate (p = 90"^) are also found to match exa_ctly with those 
given by f/tiyazaki and Silberman (1972) f'.r their H 2 -S. 

(See also 3ecs. 4,1,5 and 4,2,2). 
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RESULTS FOR AN INCLINED PLATE 

The computer code developed for the inclined plate 
solution was tested thoroughly for the hydrodynamic and 
thermal boundary layer calculations on a flat plate at 
zero incidence (j3 = 0) and for flow past a wedge. In 
these cases, of course, a similarity solution is possible 
under certain conditions (Schlichting, 1979). Results 
for the velocity and temperature profiles computed nume- 
rically agreed perfectly with the standard similarity 
solution in both cases. Following this validation of 
the computer code for the solution of momentum and energy 
equations, results were computed for the particle concen- 
tration distribution due to thermophoresis on a plate at 
zero incidence, and found to be in complete agreement 
with those of Goren (1977), These tests were done for 

various values of T , Ec and K. 

w 

Results were then computed for thermophoresis over 
a plate inclined at various angles (p = 15°, 30°, 45°, 

60° and 90°). While the Prandtl number was kept constant 
at 0.7 for most of the results, the effect of other 
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parameters such as the Eckert number, thermophore tic 
coefficient, and plate temperature was investigated. Also 
studied was the thermophoresis over an adiabatic plate. 

Most of the results are reported at five X-locations, two 
in the vicinity of the stagnation point (X = + 0.01) , and 
the other three away from the stagnation point at X = + 1.0 
and 3.0o 


4.1 Hydrodynamic Boundary Laver 

The only parameter that affects the velocity profiles 
is the pressure gradient in the momentum equation (2,28) . 
This in turn depends upon the angle of inclination, p, of 
the plate. 

4.1.1 Longitudinal Velocity Profiles 

The normalised streamwise velocity (U/U ) profiles 

o 

in the boundary layer at various angles of inclination of 
the plate (P) are shown in Figs, 4.1 to 4,3. .It may be 
noted from these figures that near the stagnation point, 

(X = 0,01) , the hydrodynamic boundary layer thickness 
increases with p while the opposite is true away from the 
stagnation point ( jxj > 1.0), There is perfect agreement 
between the present results and the standard similarity 
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solution at zero incidence (Blasius flow) denoted by * o* 
in Fig. 4.1, The velocity profiles in the negative X- 
direction are also shown in Figs. 4,2 and 4.3. These 
figures show that at a small angle of inclination, like 
P = 15°, the velocity profiles are different even when 
one moves only a small distance on either side of the 
stagnation point. The differences increase and are appre- 
ciable for larger values of p also as one moves away frcm 
the stagnation point, as expected, 

4,1,2 Transverse Velocity Profiles 

A consequence of the drastic change in the pressure 
distribution, for incidence angles other than zero, is 
more clearly reflected in the transverse velocity (V) 
distribution reported in Figs. 4.4 to 4,6. Except for the 
zero incidence angle, the transverse velocity component is 
negative over the entire characteristic length ( jxj < 0.5). 
However, at larger |X| values the transverse velocity tends 
to become positive, and as expected, for smaller angles of 
inclination first. The reason for V being negative is the 
presence of a stagnation point at non-zero angles of inci- 
dence of the plate. Near the stagnation point, the negative 
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transverse velocity increases -with the angle of incidence 
(Fig® 4,5), but away from the stagnation point at X = 1.0, 
it decreases as the angle of incidence increases (Fig, 4o6), 
The symbols fo* in Fig, 4o4 again represent the Blasius 
similarity solution thereby showing perfect agreement of the 
present results for p =0. The profiles in the negative X- 
direction are also shown in Figs. 4.5 and 4.6. A behaviour 
similar to that for the longitudinal velocity is also exhi- 
bited by the transverse velocity in the negative X-region 
of the plate. 

4ol.3 Hydrodynamic Boundary Layer Thicknesses 

The various hydrodynamic boundary layer thicknesses 
are shown in Figs. 4,7 to 4,9. The ’stars* on these figures 
correspond to the exact values for a flat plate at zero 
incidence. Clearly, the presently computed results for p = O 
are in very good agreement with the exact results. The 
integrals involved in the computation of displacement and 
momentum thicknesses are evaluated using the trapezoidal 
rule. The hydrodynamic boundary layer thickness in Fig. 4.7 
corresponds to the usual U/U = 0.99. Note that the boun- 
dary layer thickness is non-zero at X = 0 for p 0 as ex- 
pected. Also, for p = 90®, all boundary layer thicknesses 
are almost constant for jXj jC 0.5. For other values of p , 

minimum value of 6 is reached in X < 0 region, as expected. 
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The location of this minima shifts away from the stagnation 
point as p increases, except for p = 90°. For large values 
of jxj the curves for diffexent p are almost parallel except 
for p = 90°. Thus while 6ixl a constant for p = 0 

for all x>0it is a constant for non-zero p 90°) for large 
only. Also, while 6 gets larger near the stagnation 
point as p increases, it decreases with p away from the 
stagnation point except for p = 90°. Here 6 refers to any 
of the three thicknesses under consideration. 

Moie chat close xo the stagnation T'OcriL., b/L 

is of oidoi' unrly. For the isounOcAry lay^.--' cinulybis to 
be valxu, it is roguj-rc'd that 6/x uo srnali. Thus it 
appears Chat for X o, the pijsent results would oe 
valid tor Re -*“«•. At such a high Ks, the assumption that 
the slot jet be laminar would no longer hola. According 
to Hanks and Ruo (1966), the upper limit of Re for ensur- 
ing a laminar slot jet would be about 2800. Thus the 
present lesults based on the boundary layer analysis of a 
laminar slot jet are expected to hold good experiment ally 
in the region lXj>0.. 1. 

4.1.4 Flow Characteristics 

Figures 4.10 through 4.13 show the streamlines over a 
cold plate at different angles of inclination. The values 
of dimensionless stream function 'i noted on the various 
curves in these figures were computed from 

Y 

'f = / U dY . 

0 


(4.1) 
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Note that the streamline pattern is unsymmetrical about the 

Y— axis for non— zero values of p < 90° but becomes symmetrical 
for p = 90°. Note also that in these as well as in 

Figs. 4.14 to 4.16, we use the boundary layer coordinate Y 
instead of the physical coordinate y. Thus it may be di- 
fficult to visualize that the streamline pattern in Figs. 

4.10 through 4.12 corresponds to the angle p as noted. 

Figures 4.14 and 4,15 show the gas and particle velo- 
city vectors, respectively, in the vicinity of a cold plate 
inclined at p = 15°. Clearly, while the gas velocity 
vectors follow the streamlines of Fig, 4.10 (note the ratio 
of ordinate scales on Figs, 4.10 and 4.14), the aerosol 
particles fall to the cold plate in the normal direction 
due to the thermophore tic velocity component V^. Recall 
that while U and V are zero at the plate, is negative 
at the cold plate. On these figures the length of the 
velocity vectors is proportional to the magnitude of the 
velocity at the mid-point of the Vector. The ratio of the 
proportionality constant for the length of the vectors in 
Figo 4.14 to that of the vectors in Fig. 4.15 is 25. Simi- 
lar results are obtain.ed at other Inclinations of the 
plate. From the streamline plots (Figs. 4,10 through 4.13) 
the nature of isovelocity curves can be easily guessed. A 
saiiple set of isovelocity curves, U/U , is shown in Fig, 4.16 

o 

for p = 15°. The sharp turn in the curve is clearly visible 
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near the stagnation point, specially in the X < O region. 
The hydrodynamic boundary layer thickness curve for P = 15° 
in Fig. 4.7 is the sarne as the curve corresponding to 
U/U^ = 0.99 in Fig. 4.16. The isovelocity curves obtained 
for other inclinations of the plate are similar to those 
shown in Fig. 4.16. 

4.1.5 Skin Friction 

The local skin friction coefficient is calculated 
using the relation 

i = (dU/dY)y^ . (4.2) 

This parameter on both sides of the stagnation point and 
at different angles of impingement is shown in Fig. 4,17. 
At X = 0, the friction factor is zero for p 0. This 
obviously follows from the fact that for non-zero angles 
of impingement there- is no velocity gradient in the trans- 
verse direction at the stagnation point. The effect of 
potential flow velocity distribution at the plate for 
different angles, p, shown in Fig. 2.6 is very clearly 
reflected on the local skin friction coefficient. Far 
downstream from the stagnation point (X >, 3) the potential 
flow velocity at the plate U approaches the free stream 
velocity U so that the skin friction coefficient for any 

' oo 

non-zero angle of impingement tends to vary as that for 
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the flow over a flat plate at zero incidence. The latter, 
as computed presently, is also shown in Fig. 4.17, and 
matches perfectly with the exact values (Schlichting, 1979) . 
For non-zero angles of impingement the skin friction co- 
efficient, starting from a zero value at the stagnation 
point, increases and reaches a maximum value close to the 
stagnation point owing to the small boundary layer thick- 
ness and large value of in this region. The maxima 
shifts towards the stagnation point with decrease in the 
angle of impingement since the potential flow velocity 
approaches the free stream velocity more rapidly for small 
values of p (see Fig. 2,6). Particularly striking is the 
shape of the curve in X < 0 region for small values of p 
owing to the rapid increase in for such cases. The 
curve for p = 90° in Fig. 4.17 agrees perfectly with that 
corresponding to H - » in Fig. 3 of Miyazaki and 
Silberman (1972). 

4,2 Thermal Boundary Layer 

Apart from the angle of inclination of the plate, the 
thermal boundary condition on the plate, and the Prandtl 
and Eckert numbers also affect the temperature distribu- 
tion in the thermal boundary layer. Considering, _ 

air as the carrier fluid, the Prandtl number is assumed to 
be 0.7 for most of the results presented. However, the 
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effect of different Pr and Ec is also studiedo Two ther- 
mal conditions on the plate, isothermal and adiabatic 
plate, are considered. 

4.2.1 Temperature Distribution 

4.2. 1.1 Isothermal Plate Without Viscous Dissipation 

Figures 4.18 to 4.20 show the normalised dimensionless 
temperature [ 9 = ( T-1) /( T, -1 ) jdistribution for selected 
angles of inclination of the isothermal plate at different 
X-locations from the stagnation point for Pr = 0.7, Ec = O. 
The behaviour of the temperature profiles vis-a-vis the 
angle of impingement p and distance X is very similar to 
that of the velocity profiles in Figs. 4.1 to 4,3. The 
symbols 'o' in Fig. 4.18 again correspond to the exact 
solution for a flat plate at zero incidence (Schlichting , 
1979), thus validating the present results. Computations 
were also carried out for Prandtl numbers of 0.07, 0.2, 

2.0 and 7.0 and also by taking into consideration the 
viscous dissipation term with Ec = 1. Results for these 
cases were found to be quite similar to those shown in 
Figs. 4.18 to 4.20 and confirmed the fact that larger 
Prandtl numbers make the thermal boundary layer thinner. 

The thermal boundary layer thickness, 6 q , corresponding 
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to the usual value 9 = 0.01, is shown in Fig. 4.21 for 
an isothermal plate with Pr = 0.7 and Ec = 0. Similar 
results hold for other Pr and Ec values. These curves for 
6 q are similar to those for the hydrodynamic boundary 
layer thicknesses in Figs. 4.7 through 4.9. For a plate 

inclined at p = 15°, the isotherms within the thermal 
boundary layer are shown in Fig. 4.22. These isotherms 
are quite similar to the isovelocity curves shown for 
P = 15° in Fig. 4.16. The thermal boundary layer thickness 
curve for p = 15° in Fig. 4.21 is the same as the curve 
corresponding to 0 = 0.01 in Fig. 4,22. Similar isotherms 
are obtained for other angles of inclination of the plate. 

4.2. 1.2 Isothermal Plate with Viscous Dissipation 

The temperature profiles for a large value of 

Ec (= 10) are shown in Figs. 4.23 and 4.24 respectively, 

for T = 0.6 and 0 , 8 . Note that for non-zero Ec , 9 depends 
w 

upon T^. The effect of viscous dissipation is illustra- 
ted in these figures for four selected angles of inclina- 
tion of the plate between 0 to 90° at X = 0.01, 1.0 and 
3.0. It is clear from these figures that the effect of 
viscous dissipation is profound as one moves away from 
the stagnation point. A temperature maximum within the 
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boundary layer is evident at such a large value of Ec . 

Note that for < 1, the minima in 9 corresponds to the 
maxima in T. For the flat plate flow (p = C) , the dissi- 
pation of mechanical energy is too rapid for the wall to 
absorb all the heat generated and the temperature exhibits 
a maxima even at X = 0.01. Hov^ever, for larger angles of 
inclination of the plate, one needs to move sufficiently 
far from the stagnation point before the temperature pro- 
file exhibits a maxima. Also, the peak value obtained at 
larger p is smaller at any X value compared to that at 
p = 0. The effect of larger wall temperature (T^ < 1) is, 
obviously, to produce a still larger peak temperature 
value within the thermal boundary layer (see Fig. 4.24). 
However, the nature of variation of the temperature pro- 
file remains similar to that for a lower T^. 

4.2. 1.3 Adiabatic Plate with Viscous Dissipatipn 

Viscous dissipation with adiabatic wall condition 
slightly increases the thermal boundary layer thickness 
at any X-location. Figures 4.25 and 4.26 show the norma- 
lised temperature profiles for various angles at Ec — 1 
and 10, respectively. Note that the profiles are similar 
to those without viscous dissipation (Fig. 4,18) except 
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that dS/dY =0 at the plate. Different Ec value also, 

have only marginal effect on tne B profiles. However, 

it should be mentioned here that, proportional to the 
amount of viscous dissipation the plate surface will be 
getting heated up. Hence when Ec 0 , the plate will be 
quite hot ( > 1) even at a small distance away from the 
stagnation point. 

4.2.2 Nusselt Number 

The Nusselt number is calculated using the relation 
Nu = hL/k^ = (dT/dY)^^, (4.3) 

\ 

where h is the local heat transfer coefficient. Figures 
4.27 to 4.30 display the local Nusselt number (based on 
the slot jet thickness L) or the normalised temperature 
gradient in the transverse direction at the plate surface 
obtained for Pr = 0.7, 1.0, 5.0 and 10.0, and isothermal 
wall condition with Ec = 0. Here again the curves for 
g = 90° agree exactly with that for H » in Figs. 4 to 7 
of Miyazaki and Silberman (1972). For p = 90° , the 
Nusselt number variation is symmetrical about the stagna- 
tion point, and Nu at X = 0 is about twice that at 

1 /s 

X = 3.0. There is a substantial change in the Nu/Re ' 
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profiles for other values of angle The maxima in these 

profiles lies in the X < 0 region, and shifts towards the 
stagnation point as p decreases. At p = 15° for Pr = 0o7, 
the maximum value of local Nusselt number is as much 10 
times that at X = 3 (Fig. 4.27). At higher Prandtl numbers 
this ratio ls even higher (Figs. 4.2S to -.30 )o However, 
this maxima lies in the region [X|<0, 1, and is therefore 
not very meaningful for the laminar slot jet as per the 
discussion in Sec, 4.1,3. Heat transfer is enhanced 
considerably by an increase in the Prandtl 
number because a large Prandtl number implies a large ther- 
mal capacity of the fluid, enabling the fluid to carry 
more heat for the same temperature difference (t^ - t^) • 

The increase in the Nusselt number at small values of p 
is also more pronounced for large Prandtl numbers. It 
should be noted that the large change in Nusselt number 
values with respect to p occurs only in the range 
-1 < X < 1. Outside this range the Nusselt number is very 
close to that for the flow over a flat plate at zero inci- 
dence. The latter is also shown in Figs. 4.27 to 4.30 
based on present computations and matches perfectly with 
the exact values (Schlichting, 1979). 
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The mean Nusselt number over the length X of the 
plate can be computed using the relation 

Nu^ Re"^/^ dX . (4.4) 

The variation of mean Nusselt number computed over a 

length of the plate equal to the slot jet thickness on 

either side of the stagnation point is shown in Fig. 4.31 

for selected values of angle p and Eckert number. As 

expected, the mean Nusselt number increases with the 

Prandtl number. The presence of viscous dissipation with 

Ec = 1 further increases the mean Nusselt number as shown 

- 1/2 

in Fig. 4.31. Also note that on a log-log plot Nu^ Re 
varies linearly with Pr. Thus we have 

Nu Re”^/^ «= Pr^ , (4.5) 

m 

where the exponent m depends upon p and Ec , and for the 
results shown in Fig. 4.31 varies from l/3 to l/2, 

4.3 Concentration Boundary Layer 

In addition to the angle of inclination of the plate, 
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and Prandtl and Eckert numbers, the parameters that 
affect the particle concentration are the thermophore tic 
coefficient and thermal boundary condition on the wall 
(T > 1 or T < 1 or (dl/dY)^ _ = 0). While the Prandtl 
number is kept constant at 0«7 throughout most of the 
present analysis (assuming air as the carrier fluid), 
effect of all other parameters is studied. Results for 
the flat plate at zero incidence are included here for 
the sake of completeness and comparison. They match 
exactly with those obtained by Goren (1977) using the 
similarity analysis. 

4.3.1 Particle Concentration Profiles for a Cold Wall 

From Eqn. (2.30) it follows that for cold wall condi' 

tions (T < 1) the thermophore tic velocity is directed 
w 

towards the wall. In this case a non-zero concentration 
at the wall is expected from the solution of Eqn. (2.31). 
In Fig. 4,32 a number of aerosol particle concentration 
profiles are given for various angles of inclination of 
the plate; all for Pr = 0.7, = 0.25, Ec = 0 and 

K = 0.75. These profiles are drawn at X = + 0.01, + 1.0 
and 3.0, and are very similar to the U-velocity profiles 
in Figs. 4,1 to 4.3 except that unlike U, C 0 at the 
cold plate. The concentration profiles at 90® angle 
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of inclination of the plate are almost stationary with 
respect to X over jXj ^ 1. 

It is interesting to note that the wall concentration 
(at Y = O) is almost independent of the angle of inclina- 
tion of the plate and of the distance from the stagnation 
point. The reason for this could be attributed to the 
coupled variation of the temperature and concentration 
gradient at the wall for any P (see Eqn. (1.6))o For the 
conditions pertaining to Fig. 4.32 the value of wall con- 
centration varies from 0.31 to 0.32. Actual variation of 

C is discussed later, 
w 

The nature of particle concentration profiles is not 

altered by changing the wall temperature as long as 

T < 1. Fiqures 4.33 to 4.35 report the results for 
w ^ 

different T values for two representative angles of 
w 

inclination of the plate (p = 45° and 90°). Figure 4.33 

shows the profiles close to the stagnation point 

(X = i 0.01) while Figs. 4.34 and 4.35 represent the 

same away from the stagnation point (X = + 1.0 and 3.0, 

respectively). The particle concentration at the wall 

increases with T as the concentration distribution be— 

w 

comes flatter. At high T^ values (but T^ < 1) there is 
hardly any difference in the concentration profiles for 
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different angles especially at larger X values (see 
Figs. 4.3-^ and 4.35), 

The effect of Prandtl number on the heat transfer 
characteristics was presented in Section 4.2.2. In Fig. 
4.32 the particle concentration profiles were illustrated 
for a cold plate at Pr = 0.7 and with Ec =0, K =0.75 
and = 0.25. Keeping other conditions the same compu- 
tations were carried out for higher values of the Prandtl 
number. The results are presented in Figs. 4.36 to 4.38 
for Pr =1,5 and 10, respectively. Each of these figu- 
rues shows the concentration profiles for five selected 
angles (p = 0 ,15°,30'^,60°and 90°) at three X-locations 
(X = 0.01, 1 and 3). The important inference of these 
results is that the wall concentration decreases for 
larger Prandtl number fluids. This is caused by the 
larger heat transfer due to the higher temperature gra- 
dient at the wall (see Figs. 4,27 to 4.30). Also, the 
wall concentration gradient is smaller at higher Pr 
values as seen from Figs. 4,36 to 4.38 o 

The concentration boundary layer thickness correspon- 
ding to C = 0.99 is shown in Fig. 4.39, for an isothermal 
plate with Pr = 0.7, Ec = 0, K = 0.75 and T = 0.25, Simi- 

W 

lar results hold for other Pr, Ec, K and T (< 1) values. 

W 
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Clearly the curves for < 5 ^ in Fig.” 4.3"' a- Si '.ilar 
to thosG for the hydrodynamic and thermal boundary 
layer thicknesses in Figs. 4.7 and 4.21 respectively. 

Very close to the stagnation point (X — O.Ol) the 

concentration boundary layer thickness decreases as 
the angle of incidence decreases. Hov/ever, the 
situation at X 1.0 is just the opposite as 

shown in Fig. 4.39. The maximum concentration bou- 
ndary layer growth occurs for the flow past a 
flat plate at zero incidence. For the cold iso- 

thermal plate inclined at ^3 = 15*^, isoconcentration 
curves are shewn in Fig. 4.40 for Pr = 0.7, szc = 0, 

K = 0.75 and T. = 0.25. These curves are quite simi- 
lar to the isovelocity curves and isotherms shown 
in Figs. 4.16 and 4.22 respectively. The concentrat- 

ion boundary layer thickness curve for ^ = 15 in 
Fig. 4.39 is the same as the line corresponding 
+ Q 0 — go Fig. 4,40. Similar isoconcentrai-ion 

curves are obtained for other angles of inclination 

of the cold plate. 



88 


4e3»2 Particle Concentration at a Cold Wall 

In Fig. 4.41 the particle concentration at the wall, 

C, , is plotted against the thermophore tic coefficient for 
w 

Ec = 0 and Pr = 0.7. This shows clearly the effect on 
due to the presence of the stagnation point for ^ 0 . 

Any non-zero angle of inclination of the plate is represen- 
ted by a single curve for a particular irrespective of 
the X“location. Such curves are shown for = 0.25, 0.5 
and 0.9 in Fig. 4.41. Clearly, the value of at the 

cold plate increases as T increases and K decreases, 

w 

except when p = 0 and K - 0. The curves f or P 0 approach 
unity as K. approaches zero. However, for p = 0 the beha- 
viour of is quite different for small values of K. In 

this case C -* 0 as K ■* 0. When K = 0, there is no 
w 

thermophoresis and the aerosol particles move with the 
same velocity as the gas particles (V^ = 0). For p = 0, V 
is positive and hence the particles move away from the 
plate making -♦ 0. The situation is reversed for p O, 

as V is negative in which case - 1 as K 0. The wall 
concentration also varies with X in the absence of the 
stagnation point (P = O) as shown in Fig. 4.42 at X = 0.01 
and X = 1.0. The particle concentration values at the 



cola wall corresponding to Figo 4o41 and Fig. 4.42 are 
given in Table 4.1, 

The variation of wall concentration with respect to 

distance X at different angles of impingement is shown in 

Fig. 4.43 for Pr = 0.7, Ec = 0, K = 0.7 and T = 0.25. 

w 

This confirms the finding of Gokoglu and Rosner (1986) that 

the mass transfer due to thermophoresis is enhanced for 

non-zero pressure gradients in the boundary layer. The 

major variation of C is noted in 3 = 0 profile. For all 

w 

other angles the C„, value remains between 0.31 and 0.32 

w 

throughout the flow field as mentioned in the earlier 
section. However, a sudden decrease in the wall concentra- 
tion (from 0o32 to 0.31) is noted very close to the stagna- 
tion point in the negative X-region for lower angles of 
inclination of the plate. The reason for this is better 
understood by the wall concentration gradient curves plo- 
tted for identical conditions in Fig. 4,44. The concentra- 
tion gradient at the wall exhibits a sharp peak for smaller 
angles of impingement very close to the stagnation point in 
(X < 0) region. The local dimensionless concentration flux 
at the wall -(CVj)^, as defined by Homsy et al. (1981), is 
shown in Fig. 4.45. It may be noted that the nature of 
curves in Figs, 4.44 and 4.45 is quite similar to that of 



Table 4.1 Particle concentration at a cold w^ll for Pr 
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the Nusselt number in Figs. 4.27 to 4.30. 

Efforts were also made to find the suitability 
of Eqns. (1*6) and (1.7) for determining the particle 
concentration at an inclined plate. Table 4,2 gives a 
token comparison of values obtained by solving the 
particle concentration equation (2.31), and by using 
Eqns. (1.6) and (1.7) . The values seems to match quite 
well . 

4,3.3 Particle Concentration Profiles for a Hot Wall 

As mentioned earlier a hot wall (T >1) repels the 

w 

aerosol particles away from it. It is also evident from 
Eqn. (2.30) that T^ > 1 yields a positive thermophore tic 
velocity. Thus a critical layer is formed adjacent to 
a hot wall. However, the particle concentration at the 
edge of the critical layer deoends upon the product Pr.K 
being less than or greater than unity. The similarity 
analysis in Sec. 2.1 for thermophoresis over a wedge 
shows that the particle concentration approaches + «> as 
the hot plate is approached according as Pr.K ^ 1. 

Numerically the presence of the critical layer near 
a hot surface creates a problem for the solution of the 
particle concentration equation (2,31). Theoretically 
the solution of particle conservation equation should 
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yield C - + «» at the wall (Y = 0) deoer.ding upon the 
value of Pr.K. However, + «> cannot be used in computa- 
tions. A number of numerical experiments were carried 
out in this regard. Based on comparison of critical layer 
thickness with results of Goren (1977) for a flat plate at 
zero incidence, it was found that using C(X,0) = + 10 for 
Pr.K ^ 1 in the computations adequately represented the 
theoretical value (+ ») . Due to very high values of con- 
centration gradient, |dC/dY|, at the edge of the criti- 
cal layer, use of a grid adapted in the Y-direction on 
the basis of dC/dY is essential for accurate results. 
Figures 4.46 and 4.47 give a token comparison of results 
obtained for Pr.K < 1 with and without the use of adapted 
grid. Similar results were found for Pr.K > 1. Grids 
adapted on the basis of dU/dY and dl/dY were also tried 
but dc/d Y-adapted grids yielded the best comparison with 
Goren* s results for a flat plate at zero incidence. All 
the results were therefore found using dC/dY-adapted grid 
in the Y-direction. 

Figure 4.48 displays the concentration profiles at 

various locations on a hot plate for different angles of 

inclination of the plate with T = 2.0, Pr = 0,7, Ec = 0 

w 

and K = 0.75 so that Pr.K < 1. Particularly striking is 
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the sudden drop to zero of the particle concentration at 
the critical layer, which is in accord with the observed 
sharpness of the boundary of the particle free layer on 
hot objects. It is evident from Fig. 4.48, that the cri- 
tical layer thickness depends upon both the X-location and 
the angle of inclination of the plate. Similar to the 
behaviour of concentration profiles over a cold plate, the 
change in the concentration profile with respect to X is 
minimum for p = 90° and maximum for p = 0 . At higher 
wall temperatures the concentration drop is more gradual 
as also observed by Goren (1977) for p = 0. For hot wall 

conditions (T > l) » the particles are repelled from the 
w 

wall when Pr.K < 1, whereas the concentration at Y - «> is 
maintained at C = 1 by the external flow. Hence a peak 
in the concentration profile (C > 1) is obtained outside 
the critical layer and C suddenly drops to zero at the 
edge of the critical layer. The peak concentration value 
attained adjacent to the critical layer is larger for non- 
zero pressure gradient flow (p 0) than for the flow 
without the pressure gradient (p =0). However, far away 
from the stagnation point (X = 1.0 and 3.0) only a margi- 
nal difference can be noted in the maximum concentration 

value between the two cases. For T = 2.0 the maximum 

w 

concentration value obtained is about 1.18 near the 
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stagnation point. Other hot plate temperatures 
for Pr.K < 1 show similar profiles except that 

the maximum concentration, attained near the 

edge of the critical layer, increases with 
the wall temperature at all angles of 

inclination of the plate. 


Figure 4.49 illustrates the particle concentr- 


ation 

profiles 

on a heated 

wall 

for the case 

of 

Pr.K 

> 1 with 

Pr = 0.7, 

Ec = 

0, = 2.0 

and 


K = 1.5. As expected from the earlier discussion, 

the particle concentration becomes infinite as the 
wall is approached. But for this change, the 

behaviour of the concentration profiles for diffe- 
rent angles of inclination is similar to that 
discussed for Pr.K < 1. It may be noted, 

however, that Pr - 0.7 for a gas. Thus for Pr.K 
to be greater than unity, K would have to be greater 

than about 1.4. This is outside the range of 




values for K (cf. Chapter 2). Therefore the case 
of Pr.K > 1 is one of academic interest only, as long 

as the carrier fluid is a gas. 


4,3.4 Critical Laver Thickness near a Hot Wall 

Figures 4.50 and 4.51 give the thickness of the cri- 
tical layer, , near a hot wall as a function of the 
thermophore tic coefficient for three different angles of 
inclination of the plate (p =0, 15° and 90°) at the 
locations X = 1.0 and 3.0 respectively. Curves are shown 
for various values of > 1 and for Pr = 0.7 and Ec =0. 

As expected, the critical layer thickness increases with 
the wall temperature and K. Away from the stagnation 
point ( at X = 1.0 and 3.0) Y is minimum for P = 90° and 
maximum for p = 0. The variation in critical layer thick- 
ness with K and X-location (upto X = 1.0) is shown in 
Figs. 4,52 to 4.54 for three different angles, 0, 15° and 
90°, respectively. The curves are shown at X = 0.01, 0.1 
and 1.0 for three different wall temperatures (T^, = 2.0, 
1.25 and l.l) and with Ec = 0 and Pr =0.7. For a given 
K, the effect of increasing X as well as T^ on the cri- 
tical layer thickness diminishes as p increases. Tables 4.3 



Table 4.3 Critical layer thickness near a hot wall 
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Table 4.4 Critical layer thickness near a hot wall 
(p = 150 , Pr = 0o7, Ec= 0) 



Table 4,5 Critical layer thickness near a hot wall 
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to 4<,5 give the critical layer thickness, respectively, 

for p =0, 15° and 90° for five selected values of K with 

T = 2.0, 1.25 and 1.1 at X = 0.01 and 1.0. 
w 

In Figs. 4.55 and 4.56 the variation of critical 
layer thickness with respect to wall temperature is shown 
for K = 0.5 and 2.0, respectively. For the same two K 
values Figs. 4.57 and 4.58 give the variation of with 
respect to the angle of inclination of the plate. In 
these figures curves are shown for at X = 0.01 and 1.0. 
From Figs. 4.55 and 4.56 it may be noted that the critical 
layer thickness changes only a little with X over X ^ 1 
for p = 90° but for p = 0, it changes a lot. Figs. 4.57 
and 4.58 show that while Y^ increases almost linearly 

with p at X = 0.01, it decreases first sharply (for 
p < 20°) and then gradually as p increases at X = 1.0. 

4.3.5 Viscous Dissipation Effect 

The effect of viscous dissipation becomes important 
at high speeds, i.e., for Ec of the order of unity. Re- 
sults were therefore computed for both cold and adiabatic 
wall boundary conditions for Ec ranging from 1 to 10. 

> 1 boundary condition leads to results very similar to 
those for the adiabatic wall and therefore results for 
heated wall when Ec 0 are not presented. 
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4. 3. 5.1 Cold Wall 

Figures 4.59 to 4.62 illustrate the effect of vis- 
cous dissipation at four selected angles of inclination 
of the plate between 0 to 90°. These figures show the 

particle concentration profiles for Pr = 0.7, T = 0.6 

w 

and K = 0.75 at X = 0.01, 1.0 and 3.0. Near the stagna- 
tion point (X = 0.01) the profiles are almost independent 
of the Eckert number except for |3 = 0. Particularly 

striking is the drop in C as Ec is raised from 0 to 10. 

w 

As illustrated in Figs. 4, 23 and 4.24 the isothermal 
plate temperature profiles exhibit a maxima within the 
boundary layer at large values of Ec. Correspondingly, 
the temperature gradient (dl/dY) is negative in a major 
portion within the thermal boundary layer. Hence with 
non-zero viscous dissipation the thermophore tic velocity 
(directly proportional to dT/dY) becomes positive in the 
region where the temperature gradient is negative. Due 
to this the particles are blown away from the surface 
although is more or less the same at the wall for any 
value of Ec. Consequently, a significant reduction in 
particle concentration is noted near the plate along with 
a corresponding increase in values greater than unity a 
little away from the plate. The peak concentration value 


102 


obtained decreases marginally and shifts closer to the 
plate as p increases (see Figs. 4.59 to 4.62). For p = 0, 
the wall concentration remains almost constant over the 
entire plate surface for any value of Ec . However, when 
the stagnation point is present (p ^ 0 ) , the wall concen- 
tration decreases as X increases. Similar results were 
also obtained with larger cold wall temperature (T^ = 0.8). 
They are presented in Figs. 4.63 through 4.66. The curves 
for nondimensional concentration profiles in Figs. 4.59 
through 4.66 are in perfect agreement with the results of 
Goren (1977) and Gokoglu and Rosner (1985) shown in their 
Figs. 6 and 5 , respectively . It may also be noted that the 
wall concentration depends upon T only at lower values of 
Ec. Irrespective of the angle of inclination , the wall 
concentration (away from the stagnation point) is higher 

at larger T values when Ec < 3.0. 

w 

4. 3.5.2 Adiabatic Wall 

For the flow over an adiabatic inclined plate, Figs. 
4.67 to 4.70 show the particle concentration profiles for 
a series of Eckert numbers (Ec =0, 1, 3, 5 and 10) and 
angles of inclination (p = 0, 15° 45°and 90°). Here again 
the profiles are shown at X = 0.01, 1.0 and 3.0 for Pr = 0.7 
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and K = 0.75. These profiles exhibit a behaviour similar 
to that for the hot wall condition explained in Sec. 4.3.3, 
except near the wall (Y 0) where = 0 for the present 
case since (dT/dY)Y;^ = 0 leading to = 0. It can be 
observed that the effect of Ec as well as X-location on 
the concentration profile diminishes as p increases. This 
follows from the negligible variation of temperature pro- 
files over an adiabatic plate at non-zero angles of incli- 
nation for Ec 0 (Sec. 4. 2. 1,3 and Figs. 4.25 and 4.26), 
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Fig. 4.1 Longitudinal velocity profiles for flow over an inclined plate 
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Fig. 4.7 Hydrodynamic boundary layer thickness Vs X for various 
values of 
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Fig. 4.9 Momentum thickness Vs X for various values of P 
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Fig. 443 Flow field stream lines for 
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Fig. 4.14 Gas velocity vectors in the vicinity of a cold inclined plate 
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Fig. 4.17 Local skin friction coefficient over an inclined plate 









Fig, 4.18 Normalised temperature profiles over an isothermal plate 
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Fig.4.21 Thermal boundary layer thickness over an isothermal plate 
(Pr=0.7, Ec = 0) 
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Nusselt number over an isothermal plate 
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Fig. 4.32 Pariicle concentration profiles for a cold inclined plate 
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Fig. 4.39 Concentration boundary layer thickness over a cold plate 
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Fig. 4.53 Critical layer thickness as a function of K 
with Tw and X as parameters for 13= 15° 






0.8 


0.4 


Pr = 0.7 
Ec=0 


Fig. 4.57 Variation of critical layer J^ickness 
with angle of inclination at K= 0.5 



. .rUicQl layer 

=»Sh“a%le ^ 
















ates 



























Ec ® 
0-10 


168 



Fig 4.65 Particle concentration profiles for flow past a cold plate 
at different rates of viscous dissipation ( /3 =45^ = 0.8 
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CHAPTER 5 


RESULTS FOR A CIRCULAR CYLINDER 


The governing equations for thermophoresis over a 
circular cylinder in cross flow and over an inclined 
plate are identical* However, the evaluation of pre~ 
ssure gradient term in the momentum equation (2.28) is 
different in the two cases. As described in Sec. 2.4.2 
pressure gradient in the boundary layer over a circular 
cylinder in cross flow is given by Eqn. (2.42) or Eqn. 
(2.44). Thus the computer code develooed and tested 
for the inclined plate was used for the circular cylin- 
der as well after replacing the subroutine for deter- 
mining the pressure gradient. Several cases were run 
in order to cover a wide range of parameter values for 
the theoretical as well as the experimental pressure 
distribution given by Eqns. (2.42) and (2.44) respecti- 
vely. It was found that the boundary layer separates 
over the cylinder at 105.5° and 81° respectively for the 
theoretical and experimental pressure distributions. 

This is in perfect agreement with the values reported in 
literature (cf. White, 1974), 
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Results are presented here for the velocity, tempera- 
ture and concentration boundary layers with cold, hot and 
adiabatic wall boundary conditions. These profiles are 
shown versus the non-dimensional distance Y in the boun- 
dary layers at different locations over the cylinder upto 
the separation point, startinq from a location very close 
to the stagnation point. A comparison between the dis- 
tributions obtained by theoretical and experimental pre- 
ssure gradients is illustrated in each case. The Prandtl 
number is kept constant at 0.7. Also the effect of vis- 
cous dissipation on particle concentration profiles is 
presented for the cold as well as the adiabatic cylinder. 
For a hot cylinder, this effect is not shown since it is 
similar to that for the adiabatic cylinder. 

5.1 Hydrodynamic and Thermal Boundary Layers 

The hydrodynamic and thermal boundary layers over a 
circular cylinder have been investigated quite extensively 
by various authors (cf. Schlichting, 1979 and White, 1974). 
These boundary layer calculations can be carried out by 
several methods, such as the Blasius series, momentum and 
energy integral method, weighted residuals method, finite- 
difference method, etc. The results for velocity and 
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temperature boundary layers presented in this section are 
therefore not new. They are included here for the sake of 
completeness. 


5.1.1 Longitudinal Velocity Profiles 

The dimensionless longitudinal velocity (U/U^) pro- 
files at various locations ($) on the cylinder are seen 
plotted in Fig. 5.1. These are obtained for the theore- 
tical pressure distribution. The velocity profiles for 
$ > 90° possess a point of inflection and the point of 

separation is reached at ^t^ep. ~ 105.5 . Velocity dis- 
tributions obtained for the experimental pressure gradient 
at three locations on the cylinder are presented in 
Fig. 5.2, along with the corresponding profiles for theore- 
tical pressure gradient at these locations. In this case 
the separation occurs at ^sep. “ • Since the two 

pressure gradients are widely different for $ > 60 , the 

velocity profiles are also considerably different- from 
each other. 


5.1.2 Transverse Velocity Pro files 

The transverse velocity (V) profiles are shown in 
Fig. 5.3 for the theoretical pressure distribution. Just 
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like the flow over an inclined plate (Sec. 4olo2) the 
transverse velocity for flow over a cylinder is also 
negative over a substantial portion of the cylinder sur- 
face (up to 70 ) . This is due to the presence of 

the stagnation point which makes the pressure distribu- 
tion much different from that for the Blasius flow. A 
comparison of V-profiles for the theoretical and experi- 
mental pressure distributions is shovm in Fig. 5.4 at 
the same three locations ($) as for the U— velocity pro- 
file « The transverse velocity profiles for the two pre- 
ssure distributions are also quite different for f > 60*^. 

5ol,3 Normalised Temperature Profiles 

Figure 5-5 shows the normalised temperature dis- 
tribution [ 9 = (T-1)/(T -1)] obtained for the theoreti- 
cal pressure gradient in the momentum equation and with- 
out viscous dissipation for Pr = 0.7. The thermal boun- 
dary layer behaviour in this case is quite similar to 
the hydrodynamic boundary layer behavioxir. Just like the 

longitudinal velocity profiles, the temperature profiles 

o 

also have a point of inflection for # > 90 . Figure 5.6 
shows the temperature distributions at $ = I'^^bO^ and 
81° for both the theoretical and experimental pressure 
distributions. Conclusions similar to those for O/U^ 
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profiles can be made for the temperature profiles as 
well « 

The change in temperature profile due to increase 
in Ec (with constant wall temperature) is demonstrated 
in Fig. 5.7. A small valve of Ec is not sufficient to 
cause a temperature maxima within the boundary layer but, 
eventually, as Ec is increased we observe the maxima in 
the temperature profile. Note that for < 1, 9 < 0 

implies T > 1. Thus the maxima in T corresponds to the 
minima in 9 . This affects the thermophoretic deposi- 
tion of particles as we will see later. 

The normalised temperature profiles for flow past 
an adiabatic cylinder with different viscous dissipation 
rates are shown in Fig. 5.8. With the adiabatic wall 
condition the different rates of non-zero viscous dissi- 
pation have only a marginal effect on the temperature 
profiles at any location over the cylinder. 

5.2 Concentration Boundary Layer 

5.2.1 Particle Concentration Profiles 

Since the transverse velocity around the stagnation 
point is negative for flow past a circular cylinder, 
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Eqn. (3.5) or Eqn. (3.7) is used for obtaining the con- 
centration profiles depending upon the cold or hot wall 
condition respectively. The particle concentration pro- 
files are presented in this section for both wall con- 
ditions and with theoretical and experimental pressure 
distributions in the boundary layer. 

5. 2. 1.1 Cold Cylinder 

In Figs. 5.9 and 5.10 non-dimensional particle 
concentration (C) profiles are given at various loca- 
tions ($) for flow past a circular cylinder with 

K = 0.75, T. = 0.25, Ec = 0 and Pr = 0.7. Note that 
w 

like the velocity and temperature profiles, the con- 
centration profiles near the point of separation display 
a point of inflection. A non-zero wall concentration 
occurs in this case and its value is almost independent 
of the ^-location and the pressure distribution, just 
as for the cold inclined plate (cf. sec. 4.3.1), Though 
Eqns. (1.6) and (1.7) were proposed by Epstein et al. (1985) 
for natural convection flow from a vertical plate, they are 
found to hold for flovi/ past a circular cylincU.-r as well- 
Table 5.1 gives a token comparison of values obtained 
presently by solving particle concentration equation (2.31) 
with those obtained using Eqns. (1,6) and (1.7). 



180 


Table Particle concentration at a cold cylinder surface obtained from 

h<|ns. (?.31) , (1.6) and (1.7) for Pr = 0.7, Ec =0, K =0,75 




1 , * 0, 
w 

5 

T 

w 

= 0.75 



= 0.883 


( de- 
groos) 

Eqn . 

(2.31) 

bnn . 

(i.M 

Eqn 6 
(1.7) 

Eqn. 

(2.31) 

Eqn . 
(1.6) 

Eqn . 
(1.7) 

Eqn • 
(2.31) 

Eqn 0 
(1.6) 

Eqn . 

(1.7) 

0 

0.5819 

0,5821 


0.8139 

0.8146 


0.9196 

0.9247 


5 

0.5827 

0,5B2B 


0.8143 

0.8143 


0.9200 

0.9231 


20 

0.5827 

0.5827 


0.8142 

0.6145 


0.9200 

0.9231 


40 

0.5825 

0.5824 

0.5875 

0.8142 

0.8145 

0.8120 

0.9199 

0.9233 

0.9164 

60 

0.5B21 

0.5821 


0.8140 

0.8142 


0.9199 

0.9226 


80 

0.5812 

0.5814 


0.8137 

0.8142 


0.9198 

0,9230 


90 

0 . 5804 

0.5B06 


0.8133 

0o8139 


0.9197 

0.9225 


100 

0.5785 

0.57H7 


0.8126 

0.8133 


0.9194 

0,9223 



Table b.2 Local particle deposition flux over a cold 
cylinder (Pr = 0,7, Ec = 0 , K = 0o75) 


-(C v^)„ 


9 

( decrees ) 


0 

5 

20 

40 

60 

80 

90 

100 


T = 0.5 
w 

0 • 3304 
0.3278 
0.3221 
0.3038 
0.2728 
0.2267 
0.1960 
0.1511 


T = 0.75 
w 

0.1547 

0.1534 

0.1608 

0.1421 

0.1276 

0.1061 

0.0917 

0.0707 


T = 0.883 


0.0697 
0.0689 
0.0677 
0.0639 
0.0574 
0.0477 
0.04 12 
0.0318 
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Notwithstanding the almost uniform particle concentration 
at the cylinder surface, the concentration gradient at the 
cylinder surface, and the local deposition flux, -(CVj)^ , vary 
nearly oarabolically as shown in Figs. 5.11 and 5.12, respecti- 
vely, for the theoretical pressure distribution on the cylinder 
and for Pr “ 0.7l, Ec = 0, K = 0.8, and various values of T^. 
Figure 5.11 also shows the variation of the local Nusselt number 
(based on the cylinder radius) with # . Just as in the case of 
the inclined plate, the Nu curve is similar to the concentration 
gradient curve. ihe coupled variation of the temperature and 
concentration gradients at the cylinder surface causes the wall 
concentration to remain constant irrespective of uhe location 
and tiie boundary layer pressure distribution. 

In accoxoance with the results of Homsy et al. (1981'), and 
Alam and Mehrotra (1987), which are also shown in Fig. 5.12, the 
local deposition flux decreases monotonically from a maximum 
value at the stagnation point as one moves towards the point of 
separation ( f ^ 105.5° for the theoretical pressure distribut- 
ion). ‘'lote that the deposition flux at any point on the cylinder 
decreases as rises to unity. For = 1, of course, the de- 
position flux is zero everywhere on the cylinder. It may be 

pointed out that the local deposition flux obtained by Homsy 
et al, as well as by Alam and Mehrotra has to be multiplied by 
a factor of f2 before comparison with our results is possible. 

This is due to our Re being only half of theirs 


V^^hile our 
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value of the deposition flux matches exactly with that of 
Homsy et al. for = 0.5 at the stagnation point, the diff- 
erence between the two results increases as one moves away 
from the stagnation point. This is expected, however, since 
the results of Homsy et al. involve an increasing error at lar- 
ger values of $ due to the increasing inaccuracy of the tru- 
ncated Blasius series used by them for solution. In fact, 

Homsy et al. report physically impossible negative values at 
large enough values of $ due to the use of truncated Blasius 
series. The comparison in Fig. 5.12 shows an error of only 
1.66 at ^ = 20° but an error of 72.6 at ^ = 90° for 
the results of Homsy et al. Comparison with the results 
of Alam and Mehrotra is really meaningless since they do not 
provide the values of various parameters used for their 
results. The best that could be ascertained (Adam, 1987) is 
that most probably their results hold for Pr = 0.7, Ec = 0 
and K = 0.8. However, is still unknown for their 
results. Their result is shown in Fig. 5.12 only for the 
sake of completeness. Table 5.2 gives the presently com- 
puted values for the local deposition flux at different 
locations for Pr = 0.7, Ec =0, K = 0.8 and = 0.5, 

0.75 and 0.883. 
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It is interesting to note that for conditions per- 
taining to Figs. 5.9 and 5.10 (T = 0.25, K = 0.75, 

w 

Pr = 0.7 and Ec = 0) the wall concentration obtained is 
0.32 which is almost the same as that for flow over an 
inclined plate under identical condition_s Xs®6.,_Fi,g._ 4,321,.. 

Thus thermophoresis yields an almost uniform concentration 
of aerosol particles over any two-dimensional surface 
under laminar flow conditions. It is also found that the 
concentration at the wall increases almost 
linearly with the wall temperature (T^ < 1) as shown 

in Fig. 5.13. The particle concentration profiles near the 
stagnation point and just before separation (with theoreti- 
cal dP/dX) at different wall temperature values are presen- 
ted in Fig. 5.14. The point of inflection in the concen- 
tration profile at fgep “ 105.5° is observed only when 

T,., is small, 
w 

5, 2. 1.2 Heated Cylinder 

For a heated cylinder (T^ > 1) results obtained with 
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Pr = 0.7, Ec ■= 0 , K = 0«7b and T = 2.0 are illustrated 

w 

in Figs. b.l5 and 5.16 (Pr.K < 1). The presence of a 

critical layer adjacent to the cylinder surface together 

with concentration just outside this layer exceeding that 

in the free stream is demonstrated in these figures. The 

critical layer thickness increases with f as shown in 

Fig. 5.17 for the theoretical pressure distribution, and 

for T = 2.0, 1.25 and 1.1. however, the maximum particle 
w 

concentration attained adjacent to the critical layer re- 
mains more or less constant irrespective of the location 
and the pressure distribution (theoretical or experimental). 
Computations with other wall temperatures also yield similar 
results except that the critical layer thickness and the 
maximum concentration attained just outside it are different 
in each case (both increasing with wall temperature). 

5.2.2 Viscous Dissipation Effect 

Computations were also made for flow past a ci.rcular 
cylinder with non— zero Eckert numbers in the range 1 to 10, 
thereby taking into account the different rates of viscous 
dissipation in the energy equation (2o29) • Both cold wall 
and adiabatic wall boundary conditions are considered along 
with theoretical as well as experimental pressure distribution 
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in the boundary layer over the cylinder. The resulting 
concentration profiles are presented in Figs. 5.18 to 5.21. 

/ 

5 .2.2.1 Colo Cylinder 

Figure 5.18 illustrates the influence of viscous dissi- 
pation rate on the particle concentration profiles for flow 

past a cold cylinder with Pr = 0.7, K = 0.75 and T = 0.5, 

w 

The main features of Fig, 5.18 are : 

(i) No effect of Ec on the concentration profile at a 
location close to the stagnation point, and 

(ii) Away from the stagnation point ( $ = 60°) a sudden 
drop in the wall concentration as Ec is raised from 
0 to 1 . 

This drop becomes gradual as Ec is raised further. Parti- 
cularly striking is the zero concentration obtained on the 
cylinder wall at $ = 60° for Ec = 10. A maximum non- 

dimensional concentration of about 1.2 is noted in this 
profile around Y = 1.3 even though at infinity (Y ^ 3 for 
practical purposes) the concentration drops to unity. This 
implies that at larger rates of viscous dissipation the 
rapid heating makes the profiles to behave similar to 
those for flow past a heated cylinder (Sec. 5, 2. 1.2). 
Temoerature within the boundary layer increases enormously 
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esoecially at large Ec values (see Fig. 5.7). Thus at 
large cc the particles are blown away from the surface 
because ot positive value of Vj. Consequently, a drastic 
reduction in particle concentration is noted near the wall 
aloru} wi th a corresponding increase to values greater than 
unity at the end of this lean layer. The experimental 
pressure distribution in the hydrodynamic boundary layer 
gives a marginal increase in the wall concentration at 
non-zero Eckert numbers (Fig. 5.19). Consequently, the 
peak concentration achieved in this case is smaller in 
comparison to that obtained with theoretical pressure dis- 
tribution. 

5.2.2.?, Adiabatic Cylinder 

Results for the flow past an adabatic circular cylin- 
der are presented in Figs. 5.20 and 5.21 for different 
rates of viscous dissipation, the two pressure distribu- 
tions, and for Pr •- 0.7 and K = 0.75. The peak value of 
the particle concentration achieved in this case is even 
more than 1.4 for Ec = 10. The general nature of the pro- 
files is similar to that for flow past a heated cylinder 
presented in Sec. 5.2 .1 .2, except near the wall where = 0 
since dT/dvIy^ =0. A marginal increase in the critical 
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layer thickness and a sliqht decrease in the peak concen- 
tration are the main features of the particle concentration 
orotiles obtained with experimental pressure distribution 
(FiPo 5,21). The marginal changes in the concentration 
profiles with respect to Eckert number obviously follow 
from the marginal differences in the temperature distri- 
butions at non-zero Ec values (see Fig. 5.8), 
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Fig. 5. 6 Effect of pressure distribution on the normalised temperature 
distribution for flow past a cylinder 
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Fig. 5.9 Particle concentration profiles over a cold cylinder at 
different locations (Theoretical pressure distribution) 
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Fig. 5.11 Local Nussell number and concentration gradient 
at a cold cylinder surface 

























Fig. 5.19 Effect of pressure distribution on particle concentration 

profiles for flow past a cold cylinder with different viscous 
dissipaj ion rates ((j) = 6d) 








CHAPTER 6 


CONCLUSIONS 


A finite-difference solution for laminar two- 
dimensional therm op heretic flow with non-zero pressure 
gradient in the flow field is obtained. A better under- 
standing is achieved about the mechanisms that account 
for thermophoretic deposition in such situations. Two 
specific examples of flow due to the impingement of a 
slot jet on a flat plate at any angle of inclination 
(0-90 ) and cross flow over a circular cylinder are 
considered. A marching technique along with self- 
adaptive grid in the direction normal to the body surface 
is used. 

A drastic change in the behaviour of some of the 
boundary layer parameters due to the presence of the 
stagnation point are demonstrated. The presence of a 
stagnation point for non-zero angles of the plate signi- 
ficantly affects both the pressure distribution over the 
plate and the potential flow velocity at the plate. 

Also the non-zero angle causes a drastic change in the 
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transverse velocity distribution in the boundary layer. 
Near the stagnation point the hydrodynamic and thermal 
boundary layers are thinner for lower angles of impinge- 
ment but away from the stagnation point they are thinner 
for higher angles of impingement except for p = 90°. The 
local skin friction coefficient shows a large value 
close to the stagnation point at small angles of incli- 
nation of the plate. Near the stagnation point the 
local Nusselt number attains a value as high as ten 
times that far from the stagnation point for p = 15°. 
However, this ratio reduces to less than two for p = 90°. 
At zero angle of inclination of the plate the present 
procedure yields the standard and Nu values for 
Blasius flow. Also, the present results for p = 90° 
match exactly with earlier results. For other values 
of p , this is the first study. 

The particle concentration boundary layer is stu- 
died for a wide range of parameters. Its behaviour on 
a cold wall (T^ < 1) is found to be similar to that of 
the velocity and temperature profiles. The concentration 
of particles at the cold wall is found to be almost inde- 
pendent of the X-location and the value of the pressure 
gradient when T^ and K are held constant. This is 
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interesting since it implies an almost uniform thermo- 
phore tic deposition of approximately the same value in 
any external laminar flow situation as long as the wall 
temperature (less than the free stream temperature) and 
thermophore tic coefficient are constant. Moreover, the 
value of the particle concentration at the cold wall 
increases with the wall temperature. Increase in con- 
centration at the cold wall is also noted for lower 
thermophore tic coefficients. However, this is not true 
for Blasius flow for which the wall concentration first 
rises, then drops as the thermophoretic coefficient is 
reduced to zero. 

For the heated wall condition (T, , > 1) presence of 

w 

a critical layer adjacent to the wall and particle con- 
centration just outside this layer higher than that in 
the free stream are notable characteristics of the con- 
centration profiles, when the product of Prandtl number 
and thermophoretic coefficient is less than unity. When 
this product is greater than unity, a critical layer with 
very high particle concentration exists adjacent to the 
plate. The critical layer thickness is a strong function 
of the pressure gradient and the location from the stag- 
nation point. However, only a marginal change in peak 
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concentration value is noted with respect to these 
parame ters . 

The dissipation of mechanical energy in the fluid 
also has a significant effect on the particle concentra- 
tion behaviour. In the cold wall case, particle concen- 
tration at the wall drops suddenly first as Ec is increased 
from 0 to 1 and slowly thereafter reaching a zero value 
i or £c = 10 in the case of flow past a cylinder. The cold 
wall temperature profiles exhibit a maximum within the 
boundary layer at higher Ec values. Larger Ec with adia- 
batic wall conditions yields concentration profiles simi- 
lar to that for a hot wall except that C = 0 and concen- 

W 

tration gradient at the wall also vanishes. Overall, the 
effect of viscous dissipation is significant at p = 0 
(dP/dX = 0) but small at p = 90° and for the flow past a 
cylinder . 

The numerical solution procedure employed for the 
present analysis is quite general and can be used for a 
number of similar problems of engineering importance on 
thermophoresis. Variable properties and compressible flow 
can be easily considered. Analysis can also be carried 
out on similar lines for internal flow situations, for 
example, for the thermophore tic study of developing flow 



213 


in the entrance region of pipes and channels. Of course, 
in such flows, the pressure gradient is unknown apriori. 
However, methods to solve for the hydrodynamically and 
thermally developing flow in a pipe or channel (without 
thermophoresis) are well known (Hornbeck, 1973). The 
effect of radiation neglected in the present analysis can 
become a serious limitation in some cases. For example, 
radiation effects cannot be neglected for the laser modi- 
fication of the MCVD or OVD processes currently used for 
the fabrication of graded index optical fiber preforms. 
More importantly for the flows studied here, the boundary 
layer is turbulent in practical applications. The present 
method can also be used in such a case along with appro- 
priate models of turbulence. 



APPENDIX A 


T RI DI AGON AL ALGORI TWA 


The discretized form of momentum, energy and 
particle concentration equations (3.11) through (3.16) 
written for all k = 1 (1) n transverse Y~locations at 
any marching X-location (j+1) leads to a set of simul- 
taneous equations. The coefficient matrix of these 
equations have non-zero elements clustered on the three 
main diagonals. In some cases, the concentration equa- 
tion yields non-zero elements on two main diagonals only. 
All the remaining elements of the coefficient matrix are 
zero. Generally, this set of equations may be represen- 
ted in the following form® 
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The method used for solving this tridiagonal matrix 
problem is a special adaptation of the Gaussian elimina- 
tion procedure, which is commonly known as the Thomas 
algorithm (cf. Roache, 1973). To keep the size of the 
round-off errors down, Richtmyer and Morton (1967) show 
that it is sufficient that 

A}^ 1 > 0, < 0 and + Cj^j , (A.2) 

for all k from k = 1 to k = n. 

In the procedure used for solution, a series of 
matrix operations are carried out to make all the diago- 
nal elements one and the lower non-zero diagonal elements 
vanish. The resulting set of equations is of the form 
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Fpom (Ao 3) it dirsctly follows thst 


All the remaining values are calculated in sequence 
from k = n— 1 to k = 1 by back substitution. 

This method has the great advantage that all ini- 
tially zero elements of the tridiagonal band in (A.l) 
remain zero throughout all operations. As a result, 
these elements do not require storage on the digital 
computer. Thus the only storage required is for the 
tridiagonal band and the unknown and right side column 
vectors. The total storage required is therefore (5n-2) 
locations rather than (n 4-2n) required for a complete 
matrix representation. In addition, the diagonals 
Ak» and Cj^ may be represented as vectors rather than 
as elements of a matrix. This eliminates the necessity 
of using matrix notation which requires additional ope- 
rations in the system and hence additional computer 
time. The most important time saving with this method, 
however, is brought about by the fact that the number 
of necessary operations is of 0(n) , whereas for the 
usual Gaussian elimination for a filled matrix, the 
number of operations is of 0(n ). 
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A slightly improved version of the FORTRAN code of 
Hornbeck (1973) is used for solution in the present ana- 
lysis. Sets of equations with n of the order of several 

hundred are readily and accurately handled using this 
routine . 



APPENDIX B 


VARIABLE WIESH TECHNIQUE 

For computational efficiency but without sacrificing 
accuracy of the solution, a variable grid size is required 
in both X and Y-direc tions . The mesh size in the march- 
ing X-direction can be changed at will without any modi- 
fication in the discretized equations since only first 
order forward differences of error 0( A X) are involved. 
However, the real necessity for a variable mesh is evi- 
dent along the Y-direction in which variation in the 
gradients of velocity, temperature and concentration is 
quite large. Accordingly, a self-adaptive grid scheme 
(Sec, 3<,4) was employed to determine the appropriate mesh 
size along the Y-direction, The mesh size (AY) is auto- 
matically made fine in regions of large gradients and 
relatively coarse away from it by the self-adaptive grid 
scheme. This makes every AY uniquely different from 
the other. 

The modification used here in the discretized equa- 
tions due to a variable AY is essentially of the same 
nature as that described by Hornbeck (1973). A variable 
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AY presents no difficulty in the calculation of particle 
concentration (using Eqns. (3.13) through (3.16)) since 
only forward or backward first order differences of error 
0( AY) are used in these equations. However, the use of 
central differences in the Y-direction for discretizing 
the momentum and the energy equations (Eqns. (3.2) and 
(3.3) respectively in which the error is 0 (aY^)), makes 
it impossible to vary aY without some modification in 
the previously presented Eqns. (3.11) and (3.12). 


Suppose the mesh size is to be changed from a 
smaller AY 2 to a larger AY^^ at the transverse position 
corresponding to k = p, as shown in Fig. B.l(a). A 
velocity q ^ temperature q iri the case of 

energy equation) at a point AY 2 above the point k = p 
is determined by passing a parabola through the values 

interpolating for 

yielding 




+ 2 ^j+l,p4-l » 


(B.l) 


where a is the ratio of smaller to larger mesh size, 

i .e . , 

= AYj/ iVi 


a 


(B<,2) 
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Similarly, for the change of mesh size from a larger 

value ( A ) to a smaller value ( A Y 2 ) as shown in 

Fig. B.l(b) the corresoondinq relation for U is 

J+l»q 


j+1 » q 


rP3 - 2(1-0) 


p 


+ 


0-1 


U 


0+1 j+l,p+l 


(B.3) 


where o is given again by Eqn. (B.2). 


The proper value of U. , found in Eqn. (B.l) or 
Eqn. (B.3) can be used directly at every k location 
where the adjacent mesh sizes are not equal. When the 
mesh size changes from small to large, substitution of 
Eqn. (B.l) in the discretized momentum equation (3.11) 
and energy equation (3.12) gives, after rearrangement, 
respectively ; 
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The relations obtained by substituting Eqn. (B.3) 
into the momentum and energy equations, while the mesh 
size changes from a large to a small value, are identical 
to (B.4) and (B.5). These modifications in the coeffi- 
cients of the tridiagonal matrix are incorporated in the 
computer code while solving momentum and energy equations 
f or U and T respectively. 
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Fig.B.I Change of mesh size 
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